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We review attempts to estimate the influence of global cosmological expansion on local systems. 
Here 'local' is taken to mean that the sizes of the considered systems are much smaller than 
cosmologically relevant scales. For example, such influences can affect orbital motions as well 
as configurations of compact objects, like black holes. We also discuss how measurements based 
on the exchange of electromagnetic signals of distances, velocities, etc. of moving objects are 
influenced. As an application we compare orders of magnitudes of such effects with the scale 
set by the apparently anomalous acceleration of the Pioneer 10 and 11 spacecrafts, which is 
lCT 9 m/s 2 . We find no reason to believe that the latter is of cosmological origin. However, the 
general problem of gaining a qualitative and quantitative understanding of how the cosmological 
dynamics influences local systems remains challenging, with only partial clues being so far provided 
by exact solutions to the field equations of General Relativity. 
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I. INTRODUCTION 

There is by now ample evidence that our Universe is 
expanding on average. This means that on the largest 
scales one observes redshifts from structures that are in- 
terpreted as recessional motion, also called the Hubble 
flow. To first approximation, the relative velocity be- 
tween two structures grows linearly with their mutual 
distance. The constant of proportionality is the so-called 
Hubble constant, Hq, whose value is now fairly accu- 
rately measured as being clo se to 70 Km • s^ 1 • Mpc -1 , 
see e.g. ( Komatsu et all 20091) . This means that for any 
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additional mega-parsec (Mpc = 3.262 x 10 6 lightyears = 
3.086 x 10 19 km) the recessional velocity picks up an extra 
70 kilometers per second. Clearly, typical peculiar veloci- 
ties superimpose on the global Hubble flow. For galaxies 
they can be up to 1000 kilometers per second, so that 
the Hubble flow definitely dominates at distances above 
200 Mpc, i.e. above supercluster scale. I n this respect i t is 
remarkable that Hubble's classic paper (|Hubbld . ll929| ) of 
1929 plots the velocity- distance relation of extra-galactic 
nebulae only up to 2 Mpc, though it has to be added that 
in those days distances where generally underestimated, 
sometimes up to a factor of 10. 

For pedagogical purposes the global expansion is some- 
times represented by the two-dimensional balloon model, 
in which three-dimensional space corresponds to the two- 
dimensional s urface of an i n flatin g rubber balloon; see 
e.g. §27.5 in (Mis ner et all Il973l ). At each point at- 
tached to the rubber material an observer sees other 
points attached in a state of radial recessional motion, 
the faster the further they are away. This picture is used 
to stress that each point is locally (i.e. with respect to 
the local rubber material) at rest but receding from all 
other points because space in-between is itself expand- 
ing. However, this global expansion does not affect all 
structures: Local overdensities in the matter distribu- 
tion may inhibit space from expanding. In the balloon 
model of ( Misner et aLl . [l973T ) this is represented by little 
pennies being glued onto the balloon. The rubber ma- 
terial underneath the coins does not expand due to the 
stiff glue which holds it in place. The question arises 
what, in reality, are the structures corresponding to the 
coin and what dynamical mechanism provides the glue? 
It is often heard that 'bound systems' do not partici- 
pate in the global expansion, or that systems below the 
scale of galaxy clusters 'break away' from the Hubble 
flow. But what does 'bound' and 'break away' really 
mean? 1 For example, is it obvious that the Astronom- 
ical Unit is not affected by global expansion ( compare 



aims, l ike, e.g., the modeling of celestial reference 



(Krasinsky and Brumbere 



b y global expansion ( compare 
g U2004IStandishtl2004b) or ca n 



it even be, as e.g. suggested in ( Fahr and SiewertJ , |2008h . 
that our Universe is contracting on small scales while it 
expands in the large? If so, what precisely would rule the 
relation between contracting and expanding scales? 

The purpose of this paper is to review and discuss at- 
tempts that aim to make precise and answer some of these 
fundamental questions, taking due account of the dynam- 
ical laws and the kinematical framework of General Rela- 
tivity. We will emphasize the changes in kinematical rela- 
tions within time-dependent spacetime geometries, which 
seem to be widely neglected in related discussions. 

Next to being a question of fundamental inter- 
est, the raised issue also needs to be clarified 
quantitatively in connection with more practical 



frames ( Klioner and Soffel . [2005). The specific question 
of whether the global expansion has any influence 
on the local dynamics and kinematics within the 
Solar System has recently also attracted increasing 
attention in connection wit h the so- called 'Pioneer- 
Anomaly' (lAnderson et al. , 19981: Anderson et al. , 

Markwardtl 120021 : iNieto and Turvshevl 12004 : 
l2005al lbl). henceforth abbreviated by 



20021: 



Turvshev et 



PA. Here frequency-measurements in Doppler tracking 
arc translated into standard kinematical quantities, 
like velocity and acceleration. The result shows an 
anomalous acceleration of the Pioneer satell ites directed 
towa rds the center of the Solar System. In ( Markwardtl 
120021) the magnitude of this acceleration is reported 
to be a = 8.6 ± 1.34 x 10~ 10 m- s~ 2 . Note that such 
an apparently small acceleration amounts to variations 
in spatial localization of nearly 500 kilometers after 
10 years. It so happens that the magnitude of this 
acceleration is very close to the product of the current 
value of the Hubble constant, Hq, and the velocity of 
light in vacuum: 

H ck (70 km -s- 1 - Mpc" 1 ) (3 x lO^m-s" 1 ) ^ 



7 x 10 



-10 . 



Whether this 'almost coincidence' of numbers does in- 
deed have any deeper significance can and should only 
be decided on the basis of reliable estimates within the 
dynamical framework of General Relativity. There al- 
ready exist various speculations and claims in the litera- 
ture that try to attribute the PA to either simp le kine- 
matical (e.g. ( Rosales and Sanchez-Gqm"e3 , ll998l) ) or dy- 
namical (e.g. (jFahr and Siewer tl 120081) 1 effects of a time 
varying background geometry, though none of them does 
justice to the requirements posed by General Relativity. 2 
This is clearly a very difficult task: There is very little 
analytical knowledge of how to model in terms of exact 
solutions, or at least in terms of controlled approxima- 
tions to exact solutions, the hierarchy of mutually embed- 
ded systems: Solar System — > Galaxy — > Local Group — > 
Cluster — » Supercluster — > Standard-Cosmological Solu- 
tion. Usually we expect each such system to define a 
typical length scale beyon d which we may consider it as 
quasi isolated (jCoxl . [20071 ). But, clearly, whether this is 
a valid assumption or not can only be decided on the 
basis of a self-consistent dynamical consideration. In our 
context all this suggests to first study the influence of cos- 
mic expansion on the most simple systems immersed in 
an otherwise homogeneous cosmological background. We 
will see that this already poses a number of non-trivial 
analytical as well as conceptual problems. 

In this article we will derive upper bounds for various 
effects of global expansion on local systems in the con- 



1 For a rec ent discussion on th e meaning of 'joining the Hubble 
flow' see dBarnes et ^.1,120061) . 



2 The reader will soon find out that we disagree with all such 
claims. 
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text of such simple models. The idea here is that the 
upper bounds so derived will a fortiori be upper bounds 
in more realistic models, since a further embedding of 
the system we consider into a higher structure of local 
overdensities will further suppress the influence of cosmo- 
logical expansion. This is evidently true in situations in 
which the spherically symmetric Einstein-Straus model 
applies, but can also be argued for as a result of taking 
into account small-scale anisotropies in the matter distri- 
bution, as has been done fr om first-order perturba t ions o f 
the Newtonian equations ( Dominguez and Gaitd . l200lh . 
We conclude form this that if we find the relevant upper 
bounds to be outside current experimental reach, this will 
maintain to be the case in more realistic contexts. 



which allegedly shows the absence of any relevant dy- 
namical effect of global expansion. Its original form only 
involved the dynamical action principle together with 
some simple scaling argument. Since this reference is 
one of the most frequently cited in this field, and since 
the simplicity of the argument (which hardly involves 
any real analysis) is definitely deceptive, we give an in- 
dependent treatment that makes no use of any hypo- 
thetical scaling rules for physical quantities other than 
spati al lengths and times . Our treatment, which fol- 
lows ( Carrera and Giuhnl . |2005j), also reveals that the 
original argument by Dicke & Peebles is insufficient to 
discuss leading order effects of cosmological expansion. 
It is therefore also i neffec tive in its attempt to contra- 
dict (|Pachneri . [l963lll964D . 



II. STRATEGIC OUTLINE AND RESULTS 



A. Improved Newtonian equations 

The strategies that so far have been followed are 
twofold: Either one studies modified Newtonian or spe- 
cial relativistic equations of motions for two point- 
particles with a force of mutual attraction (gravita- 
tional or electromagnetic). The modifications are de- 
rived from putting the system into a fixed standard- 
cosmological background (usually spatially flat) without 
back-reactions being taken into account. We shall discuss 
this appro ach in Sectio ns IIIII and [IVJ Our discussion, 
based on ( Carrera and Giulinil. 12005), compleme nts the 
perturbative analysis in (|Cooperstock et al. , 1998) which 
misses all orbits which are unstable under cosmological 
expansion (which do exist). In this respect we follow a 
ver y similar stra tegy as, e.g., in the more recent papers 
by ( Pried. l2005f) ( the basic id ea of which goes back at 
least to (jP^chned, Il963l Il964) 1 and also (jAdkins et all 
2007), though we think that there arc also useful differ- 
ences. We also supply quantitative estimates and clarify 
that the improved Newtonian equations of motion are 
written in terms of the right coordinates (non-rotating 
and metrically normalized). The purpose of this model 
is to develop a good physical intuition for the qualitative 
as well as quantitative features of any dynamical effects 
involved. 

Eventually the Newtonian model just mentioned has 
to be understood as a limiting case of a genuinely rel- 
ativistic treatment. For the gravitational case this is 
done in Section [V] (an alternative and more geometric 
derivation is given in Section IVI.Bp , where we employ 
the McVittie metric to model a spherically symmetric 
mass embedded in a spatially flat Friedmann-Lemaitre- 
Robertson- Walker (FLFW) universe. The geodesic equa- 
tion is then, in a suitable limit, shown to lead to the 
i mproved Newtonian mode l discussed above (see also 
(ICarrera and Giulinil . I2005h 1. The same holds for the 
electromagnetic case, as we show in Section IIVI There 
we take a slight detour to al so reconsider a cla s sic ar - 
gument by Dicke & Peebles ( Dicke and Peebles! . 1 1964 ). 



B. Exact solutions 

The other approach consists of finding exact solutions 
to Einstein's field equations for an inhomogeneous situa- 
tion that, in the most simple case, models a single, quasi- 
localized, non-rotating, electrically neutral inhomogene- 
ity within a FLFW universe. Using this inhomogeneous 
solution as background one can then study the motion 
of test particles (following geodesies in the background 
geometry) and, in particular, the influence of expansion 
on this motion. 

This approach can be subdivided into two strategies. 
The first tries to literally construct a new exact solution 
out of two known ones, so that the new solution con- 
tains a connected piece from each of the two old ones 
as isometric submanifolds. These we refer to as matched 
solutions. This is relaxed in the second, more general 
strategy, where the new solution is merely required to 
somehow approximate the relevant part of each of the two 
old solutions in some region. These we refer to as melted 
solutions. Needless to say that melted solutions offer a 
much greater variety for construction than matched ones. 
However, it is also true that often not much is known 
about the proper physical interpretation of the former. 
In this respect the matching solutions usually provide a 
much clearer picture. 

According to the above requirements, in both cases we 
shall restrict attention to spherically symmetric space- 
times which, loosely speaking, approximate a FLRW 
solution of standard cosmology for 'large radii' and a 
non-charged, non-rotating compact object characterized 
by the exterior Schwarzschild solution for 'small radii'. 
(Clearly there must be some characteristic radius in 
terms of which 'large' and 'small' radii are defined.) Also, 
one often restricts attention to the spatially flat FLRW 
models for simplicity, which also seems justified in view 
of current cosmological data which are compatible with 
spatial flatness. 
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1. Matched solutions 

A first approach to the matching idea was in i tiated 
by Einstein and Straus ( Einstein and Strau3 . ll945lll946h 
in 1945 and l ater worked out i n more analytical detail 
by Schiicking (jSchuckinel Il954[ ). Here the matched so- 
lution is really such that for radii smaller than a cer- 
tain matching radius, R v (henceforth called the vac- 
uole or Schiicking radius), it is exactly given by the 
Schwarzschild solution (exterior for a black hole, exterior 
plus interior for a star) and for radii above this radius 
it is exactly given by a FLRW universe for dust matter 
without cosmological constant (this can be generalized, 
see below). The radius R v is a function of the central 
gravitational mass M and the cosmological mass-density 
g, through the latter of which it also depends on the cos- 
mological time t. It is determined by 

4tt 



— RZ ■ q = M . 



(2) 



This formula holds for flat as well as curved FLRW mod- 
els if 'radius' is taken to mean 'areal radius', the defini- 
tion of which is that a two-sphere of areal radius R has 
a proper surface area of 4ttR 2 . In flat space the areal 
radius coincides with the proper radius (the geodesic dis- 
tance between the center and any point on the sphere), so 
that 4^ i? 3 is just the proper volume inside the sphere of 
radius R v (cf. Section IV. A[) . However, in backgrounds of 
positive (negative) curvature this expression is smaller 
(larger) than the proper volume (the proper volume 
grows faster (slower) with areal radius) and hence, for 
given g, the left-hand side of is also smaller (larger) 
than the proper mass of the dust contained within a 
sphere of areal radius R v . 

Here we recall that the gravitational mass of a lump of 
matter is not just proportional to the amount of matter 
(baryons) in that region. For example, the kinetic energy 
as well as the gravitational binding energy also contribute 
to the gravitational mass. This is expressed in formula 
(|D47[) of Appendix ID. 41 where further explanations will 
be provided. As is well known, the mathematical char- 
acterization of appropriate notions of quasi-local gravita- 
tional mass that would apply to general spacetimes is a 
notoriously difficult pro blem to which v arious attempts 
for solutions exist; see (ISzabadosL l200i for the current 
status. However, in the spherically-symmetric case, to 
which we restrict attention, the so-called Misner-Sharp 
energy gives a satisfying and convenient concept of ac- 
tive gravitational mass. Its definition will be given in 
Section IV. Al and more details, including its equality in 
value to the Hawking mass, are discussed in the Appen- 
dices E3 and El 

The original construction by Einstein and Straus and 
its analytical completion by Schiicking were quite com- 
plicated. We will give a much simpler and conceptually 
clearer description in Section [V. CI using a suitable refor- 
mulation of the condition for the matching of solutions. 
However, it is not hard to gain some intuitive under- 
standing for the matching construction and the value of 



R v as defined by ([2]). Let us for the moment restrict to 
the spatially flat case and consider the homogeneous and 
isotropic dust-filled universe at some moment of time t. 
The dust within a 3-ball of proper radius R v represents an 
amount of matter of total mass M as given by ([2]). Now 
compress this amount of matter in a spherically sym- 
metric fashion until it becomes a compact star or a black 
hole. In Newtonian gravity the gravitational field outside 
a spherically symmetric mass distribution only depends 
on the total mass and not on its radial density distri- 
bution. This is also true in General Relativity, which 
is essentially the content of Birkhoff's theorem. 3 Hence 
the above compression preserves equilibrium (albeit an 
unstable one, see below) for the dust particles just out- 
side the boundary-sphere of radius R v . For radii smaller 
than R v we have the Schwarzschild solution (which is 
the unique non-trivial spherically symmetric vacuum so- 
lution according to Birkhoff's theorem) which therefore 
matches to the FLRW solution for R > R v at the bound- 
ary R — R v where the matter density is discontinuous. 
The spatial two-sphere R = R v is comoving with the 
Hubble flow, meaning that its proper surface area grows 
in case of expansion. Finally, in case of constant positive 
(negative) spatial curvature, ((2|) tells us that the matched 
Schwarzschild solution has a smaller (larger) mass than 
the mass that the amount of dust represents within the 
ball of areal radius R v within the FLRW universe. 

The Einstein-Straus model can be generalized in sev- 
eral ways. Instead of cutting out one ball, one can 
cut several non-overlapping ones and fill in the in- 
teriors with Schwarzschild geometries of appropriate 
masses. For obvious reasons these are sometimes re- 
ferred to as 'Swiss-Cheese models'. These, in turn, 
can be generaliz ed to the cases of non -vanishing cosmo- 
logical co nstant dBalbinot et all Il988l ) or non- vanishing 
pressure (|Bona and Stelal ." 19871 ). Finally, the Einstein- 
Straus model can be generalized to spherically symmetric 
but inhomogeneous Le maitre-Tolman- Bondi (LTB) cos- 
mological backgrounds ( BonnorL l2000f) . 

Since for the Einstein-Straus model the geometry 
within R < R v is exactly Schwarzschild (for vanishing 
cosmological constant) or Schwarzschild-de Sitter space- 
time (for non- vanishing cosmological constant), it is clear 
that any dynamical system situated in this background 
geometry (no back reaction) only detects that part of the 
cosmic expansion that is due to a non-vanishing cosmo- 
logical constant. In particular, for vanishing cosmological 
constant, the cosmic expansion that goes on outside the 
expanding vacuole R = R v is not felt from within. Hence 
global expansion due to ordinary (localizable) matter 
can, in principle, be completely inhibited by local inho- 
mogeneities. 

There are, however, several severe problems concern- 



3 An elegant proof of Birkhoff's theorem will appear as a by- 
product from our considerations in Appendix ID. 31 
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ing the Einstein-Straus approach. First of all, it cannot 
provide a realistic model for the environment of small 
structures in our Universe, 'small' meaning below the 
scales of galaxy clusters or superclusters. To see this, 
apply ([2]) to a spatially flat universe whose background 
matter density g is given by the critical density 



£>crit 



3Hl 
8nG ' 



where G is Newton's constant. Then ^ gives 



R v = {Rs R 2 H ) 1/3 ^ (-£-) 4001y, 



M 

M 



1/3 



where 



Rs 
Rh 



2GM 

c 

#0 " 



M 



3 km . 



4Gpc fa 1.3 x 10 23 km, 



(3) 

(4) 

(5) 
(0) 



are the Schwarzschild radius for the mass M and the 
Hubble radius, respectively. M© — 2 x 10 30 kg is the 
solar mass. 

For a single solar mass this gives a vacuole radius of 
almost 400 lightyears, which is almost two orders of mag- 
nitude larger than the average distance of stars in our 
Galaxy. Therefore, the Swiss-Cheese model cannot ap- 
ply at the scale of stars in galaxies. This changes as 
one goes to larger scales. For example, the Virgo cluster 
is estim ated to have a mass of approximately 10 15 solar 
masses (jFouque et all l200ll ) 4 . which makes its vacuole 
radius 10 5 times larger than that for a single solar mass, 
so that it is approximately given by 10 Mpc. This is just 
a little smaller than the average distance of groups and 
clusters of galaxies within the Virgo supercluster. Hence 
the Einstein-Straus approach might well give viable mod- 
els above cluster scales. Similar conclusion can be drawn 
for th e vacuole construction in LTB spacetimes ( Bonnorl . 
I2000T) : There it is argued that the vacuole might be as 
big as the Local Group. 

The Einstein-Straus solution (as well as its general- 
ization for LTB spacetimes given by Bonnor) may also 
be criticized on theoretical grounds. An obvious one 
is its dynamical instability: slight perturbations of the 
matching radius to larger radii will let it increase with- 
out bound, slight perturbations to smaller ra dii will let 
it col lapse. This c an be proven fo rmally (e.g. ( Krasiriskil 
1 1998(1 . Ch. 3 and (|Bonnoii 120000 ) but it is also rather 
obvious, since R v is defined by the equal and opposite 
gravitational pull of the central mass on one side and 
the cosmological masses on the other. Both pulls in- 
crease as one moves towards their side, so that the equi- 
librium position must correspond to a local maximum 



of the gravitational potential. Another criticism of the 
Einstein-Straus solution concerns the severe restrictions 
under which it may be genera lized to non spherically- 
symmetric situations; see e.g. ( Mena et all 120021 120031 
l2005t ISenovilla and VeraL Il997h . 



2. Melted solutions 

The above discussion shows that the Einstein-Straus 
approach does not give us useful information regarding 
the dynamical impact of cosmic expansion on structures 
well below the scales of galaxy clusters. For this rea- 
son other exact solutions are sought. In this respect we 
wish to remind the reader on the following general as- 
pect: In physics we are hardly ever in the position to 
mathematically rigorously model physically realistic sce- 
narios. Usually we are at best either able to provide 
approximate solutions for realistic models or exact so- 
lutions for approximate models, and in most cases ap- 
proximations are made on both sides. The art of physics 
then precisely consists in finding the right mixture in each 
given case. However, in this process our intuition usually 
strongly rests on the existence of at least some 'nearby' 
exact solutions. Accordingly, one seeks exact solutions 
in General Relativity that, with some degree of physi- 
cal approximation, model a spherically symmetric body 
immersed in an expanding universe. However, it is not 
as easy as one might think at first to characterize 'body' 
and 'immersed'. 5 Clearly it is associated with some inho- 
mogeneity in form of a spatial region with an overdense 
matter distribution, as compared to that of the approx- 
imately homogeneous distribution far out. But a body 
should also be quasi-isolated in order to be distinguish- 
able form a mere local density fluctuation with smooth 
transition. Typical exact solutions that models the latter 
are the LTB solutions, in which matter is represented by 
pressureless dust that freely falls into the local overdense 
inhomogeneity. In some sense, these form the other ex- 
treme to the Einstein-Straus solutions in that they make 
the transition as smooth and mild as one wishes. Here 
we shall be interested in models that somewhat lie in- 
between these extremes. 

An attempt to combine an interior Schwarzschild so- 
lution (representing a sta r) and a fl at FLRW universe 
was made by Gautreau ( Gautreau] . Il984f ). Here the 
matter model consists of two components, a perfect 
fluid with pressure and equation of state p = p(g) out- 
side the star, and the superposition of this with the 
star's dust-matter inside the star. However, Gautreau 
also made the assumption that the matter outside the 



4 Their considerations are based on a LTB model for the cluster. 



In a linear theory, the 'simultaneous presence' of two structures, 
like a local inhomogeneity in an 'otherwise' homogeneous back- 
ground, naturally corresponds to the mathematical operation of 
addition of the corresponding individual solutions. In a non- 
linear theory, however, no such simple recipe exists. 
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star moves on radially infalling geodesies, which is 
only consistent if the pressure outside is spatially con- 
stant. T hus one is reduced to exact F LRW outside 
the star (jvan den Bergh and Wilsl 1 1984 or the LTB 
mode l. (Further remarks may be found in I Krasihskll 
Il998l ). e.g. p. 113 and 165.) Other solutions, model- 
ing a black hole in a cosmological spacetime, have been 
given in the literature. However, these solutions model 



objec t s which are ei t her ro tating (jRamachandra et 



2003; IVaidval . 



2004), or both 



197 1 119841). charged (IGao and Zhane 



Patel and Trivedil . ll982l) . Surveys on the 



subje ct of c osmologic a l blac k holes are ( Vishveshwaral 
l2000h and (|McClurd . 120061). Further int e restin g so- 



lut ions are given in ( Raiesh Navak et al. j_ 200 ll) an d 



( R^^ 

in (jFaraoni and Jacquesl 12007^ Sultana and Dverl . l2005h . 



The solutions proposed in the latter two wor ks can be 
seen as generalizations of McVittie's model ( McVittiel . 
1933), which we extensively discuss in Section IV. Dl A 
crucial feature of these solutions is, however, that the 
strength of the inhomogeneity 6 varies in time, whereas 
for the McVittie model it remains constant. These 
solutions are of in terest in their own right (for a de- 
tailed analysis see ( Carrera and Giulinil . [20091). but our 
goal here is to focus on the effects due to cosmologi- 
cal expansion and not on the effects due to a chang- 
ing strength of the central inhom ogeneity. The solu- 
tion proposed in the former work (jRaiesh Navak et all 
2001) is the melting of a Schwarzschild spacetime in 
an Einstein's static universe. This is a purely static 
sol ution whose properties and geodesies where studied 
in ( Ramachandra and Vishveshwaral |2002[ ) . For our pur- 
poses, however, this spacetime is not interesting since it 
is asymptotically an Einstein universe, and hence not in 
agreement with the present picture of our Universe at 
large scales. 



For these reasons in Section IV. Dl we shall pay spe- 
cial attention to the McVittie model. This contains a 
distinguished central object in the sense that the mass 
within a sphere centered at the inhomogeneity splits into 
a piece that comes from the continuously distributed cos- 
mological fluid (with pressure) and a constant piece that 
does not depend on the radius of the enclosing sphere; 
see our Eq. (f8"2"| . Moreover, the latter piece is also con- 
stant in time, meaning that the strength of the central 
inhomogeneity remains constant. By the way, McVit- 
tie's solutions contain the Schwarzschild-de Sitter one as 
a special case, which was recently used in the litera- 
ture to estimate t he effects of cosmological expansion 
on local systems (Hackma nn and Lammerzahll . |2008c; 



Kagramanova et al. 



20061 ). In Section lYTXlT we show 



that in a suitable weak-field and slow-motion approxima- 
tion the geodesic equation in McVittie spacetime reduces 
to the improved Newtonian equations discussed earlier. 



An alternative and more geometric derivation of the im- 
proved Newtonian equations for the McVittie case is pre- 
sented in Section IVLBl (see Eq. JUS])). 



C. Kinematical effects 

1. Timing and distances 

Neither the improved Newtonian model nor other gen- 
eral dynamical arguments make any statement about 
possible kinematical effects, i.e. effects in connection with 
measurements of spatial distances and time durations in a 
cosmological environment whose geometry changes with 
time. This is an important issue if one wants to perform 
the tracking of a spacecraft, that is a 'mapping out' of its 
trajectory, which basically means to determine its simul- 
taneous spatial distance to the observer at given observer 
times. But we know from General Relativity that the 
concepts of 'simultaneity' and 'spatial distance' are not 
uniquely defined. This fact needs to be taken due care 
of when analytical expressions for trajectories, e.g. solu- 
tions to the equations of motion in some arbitrarily cho- 
sen coordinate system, are compared with experimental 
findings. In those situations it is likely that different kine- 
matical notions of simultaneity and distance are involved 
which need to be properly transformed into each other 
before being compared. For example, these transforma- 
tions can result in additional acceleration terms involving 
the product fTJ. Accordingly, there were claims in the 
literature that these kinematical effects could a c count 
for the PA; see e.g. dNieto et all 120051; iNottald. 12003 : 
Pallel 120051: iRanadaL 120051: iRosales and Sanchez-Gomea . 
ll998HRosalesl 2002) and also statements to the contrary 



( Lammerzahl et all 20061). In Section IVI.AI following 



( Carrera and Giulinil . I2005D . we will confirm the exis- 
tence of kinematical acceleration terms proportional to 
Hqc, but they are suppressed with additional powers of 
[3 — v/c, which renders them irrelevant as far as the PA 
is concerned. 



2. Doppler Tracking 



In Section lV.D.ll we will identify the strength of the inhomogene- 
ity with the Weyl part of the Misner-Sharp energy. 



The discussion in Section IVI.BI is based on 

( Carrera and Giulinil . I2006D . We explain in some de- 
tail the geometric theory for setting up the kinematical 
framework in which Doppler tracking should be discussed 
in order to properly speak of relative velocities and ac- 
celerations. This is a non-trivial issue which is, in our 
opinion, not properly appreciated in the literature on this 
subject (rela ted general discussions are (|Bini et a/.l . ll995l: 
lBolosll2007l) ). Using this setting, we show how to derive 
an exact Doppler-tracking formula for a flat FLRW uni- 
verse. This we use to give reliable upper bounds for kine- 
matical effects caused by cosmic expansion. We also dis- 
cuss generalizations to McVittie spacetime. Even though 
such effects exist, they again turn out to be irrelevant for 
the PA. 
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III. NEWTONIAN APPROACH 

In order to gain intuition we consider a simple bounded 
system, say an atom or a planetary system, immersed 
in an expanding cosmos. We ask for the effects of this 
expansion on our local system. Does our system expand 
with the cosmos? Does it expand only partially? Or does 
it not expand at all? Here we shall not be concerned with 
the far more complex problem of how stable large-scale 
structures may emerge from unstable local gravitational 
dynamics in an expanding universe. Thi s has been dis- 
cussed in (Bu chert and Domingued . l2005f ) and references 
therein. 



A. Restricted two-body problem in an expanding universe 

We consider the dynamical problem of two bodies at- 
tracting each other via a force with 1 /R 2 fall-off. For sim- 
plicity we may think of one mass as being much smaller 
than the other one, though this is really inessential. One 
may think of two galaxies, a star and a planet, a planet 
and a spacecraft, or a (classical) atom given by an elec- 
tron orbiting around a proton. The system is placed into 
an isotropically expanding ambient universe. We wish to 
know the leading order influence of the ambient expan- 
sion onto the relative two-body dynamics. 

To leading order, the global expansion is described by 
the simple linear Hubble law, R = HR, which states that 
the relative radial velocity of two comoving objects at a 
mutual distance R grows proportional to that distance. 
More precisely, the term 'distance' is here understood 
as the geodesic distance in the spacetime hypersurface 
of constant cosmological time t between its two intersec- 
tion points with the two worldlines of the objects consid- 
ered. H denotes the Hubble parameter, which generally 
depends on t but not on space. It is given in terms of the 
scale parameter, a(t), via H = a/ a. 

Taking into account H = (a/a) — H 2 , the acceleration 
that results from the Hubble law is simply given by 



R\ 



cosm.acc. 



HR + HR 



where 



aa 

n ^ 



-R 



a -H- 2 
a 



-qH 2 R, (7) 



(8) 



is the dimensionless deceleration parameter. To get a 
feeling for the magnitude, we remark that for the cur- 
rent best-estimates for the parameters H and q, Hq ~ 
70 Km • s _1 • Mpc~ x and qo w —0.6 respectively, we get 
a/ a w 3 x 10~ 36 s~ 2 , which even at Pluto's distance of 
40 AU merely amounts to a tiny outward pointing accel- 
eration of 2 x 10~ 23 m • s -1 . 

Now note that, in the sense of General Relativity, a 
body that is comoving with the cosmological expansion 
is moving on an inertial trajectory, i.e. it is force free. 
On the other hand, according to Newton, a dynamical 



force is, by definition, the cause for deviations from iner- 
tial motion. In the present context this would mean that 
dynamical forces are the causes for deviations from the 
motions described by ([7]), which suggests that in New- 
ton's law, mx = F, we should make the replacement 



R^ R- (a/a)R 



(9) 



in order to apply to the (sufficiently slow) motion of in- 
teracting point masses in an expanding universe. Note 
that this also applies to gravitational interactions in a 
Newtonian approximation in which gravity is considered 
to be a force in the above sense. 

As we will see, the replacement © can be justified rig- 
orously in a variety of contexts, like for gravitationally 
bound systems, using the equation of geodesic deviation 
in General Relativity. Whenever we attempt to justify 
the replacement © we must not forget that the Newto- 
nian equations of motion (without Coriolis and centrifu- 
gal type 'forces' in them) necessarily refer to preferred 
systems of coordinates which are 1) locally non rotat- 
ing, 2) whose origin is freely falling, and 3) in which the 
coordinate values directly refer to (local) inertial time 
(time coordinate) and spatial geodesic distance (space 
coordinates), as measured by comoving clocks and rods. 
This is achieved by using so -called Fermi norma l coor- 
dinates (see, e.g., §13.6 of ( Misner et all Il973f) ) in a 
neighborhood of a geodesic worldline — e.g. that of the 
Sun or the proton. This is also the approach followed 
in ( Cooperstock et all . Il998l ). Note that a Fermi sys- 
tem of coordinates can be defined for worldlines of arbi- 
trary acceleration and correspond to locally non-rotating 
frames, which may physically be realized by a system of 
at least two non collinear gyros in torque-free suspen- 
sions taken along the worldline. Along a geodesic, that is 
a worldline of zero acceleration, this system corresponds 
to a local inertial observer and is called Fermi normal. 
The equation of geodesic deviation in these coordinates 
now gives the variation of the spatial geodesic distance 
to a neighboring geodesically moving object, e.g. a planet 
or spacecraft. It reads 7 



d 2 x k 
~aW 2 



R k , 



M0 J 



0. 



(10) 



Here the x k are the spatial non-rotating normal co- 
ordinates whose values directly refer to the proper 
spati al distance. In these coordinates we further 
have ( Cooperstock et all 11998) 



R> 



-5fd/c 



(11) 



on the worldline of the first observer, where the overdot 
refers to differentiation with respect to the cosmological 



7 By construction of the coordinates, the Christoffel symbols 
vanish along the worldline of the first observer. Since this world- 
line is geodesic, Fermi- Walker transportation just reduces to par- 
allel transportation. 
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time, which reduces to the proper time along the ob- 
server's worldlinc. 

Equations (| 1011 1 [) simply state that in Fermi normal 
coordinates around one inertial observer another nearby 
inertial observer is radially accelerating away at a magni- 
tude (d/a)R, just as envisaged before. In linear approxi- 
mation, this acceleration has to be added to that result- 
ing from the other metric perturbation that is caused by 
the mass at the position of the first observer. As a re- 
sult, in the case of purely gravitational interaction, wc 
obtain the equation of motion of a test particle (whose 
metric perturbation we neglect) in the gravitational field 
of a heavier object whose metric perturbation away from 
the FLRW cosmological background we approximate to 
linear order. In the case of charged objects, we neglect 
the metric perturbations caused by the masses of both 
charges as well as their electromagnetic field, and simply 
take into account their mutual electromagnetic interac- 
tion. 

Neglecting large velocity effects (i.e. terms quadratic 
or higher order in v /c) we can now write down the equa- 
tion of motion for the familiar two-body problem. After 
specification of a scale function a(t), we get two ODEs 
for the variables (R,(p), which describe the position 8 of 
the orbiting body with respect to the central one: 



R 2 ip = L . 



(12a) 
(12b) 



These are the (a/a)-improved Newtonian equations of 
motion for the two-body problem, where L represents 
the (conserved) angular momentum of the planet (or elec- 
tron) per unit mass and C the strength of the attractive 
force. In the gravitational case C = GM, where M is 
the mass of the central body, and in the electromagnetic 
case, for the electron-proton system, C = e 2 /m (Gaus- 
sian unit), where e and m are the electron's charge and 
mass, respectively. In Sections [V] and IIVI we will show 
how to obtain (|12p in appropriate limits from the full 
general relativistic treatments. 

We now wish to study the effect the a term has on the 
unperturbed Kepler orbits. We start with the obvious 
remark that this term results from the acceleration and 
not just the expansion of the universe. 

Next we point out that in the concrete physical cases 
of interest, the time dependence of this term is negli- 
gible to a very good approximation. Indeed, putting 
/ := a/a, the relative time variation of the coefficient 
of R in ([7]) is ///. For an exponential scale function 
a(t) oc exp(Ai) (vacuum-energy-dominated universe) this 
vanishes, and for a power law a(t) cx t x (for exam- 
ple matter-, or radiation-dominated universes) this is 



Recall that 'position' refers to Fermi normal coordinates, i.e. R 
is the radial geodesic distance to the observer at R = 0. 



—2H/\, and hence of the order of the inverse age of the 
universe. If we consider a planet in the Solar System, 
the relevant time scale of the problem is the period of its 
orbit around the Sun. The relative error in the distur- 
bance, when treating the factor a/a as constant during 
an orbit, is hence smaller than 10 -9 . For atoms it is much 
smaller, of course. In principle, a time varying a/ a causes 
change s in the semi-major axis a nd eccentricity of Kepler 
orbits ( Sereno and Jetzerl |2007| ). But here we shall ne- 
glect the time-dependence of ([7]) and set a/ a equal to a 
constant A. Because of © we have A := —q H 2 . Then 
(j!2a[) can be immediately integrated: 



l -R 2 + U{R) = E, 

where the effective potential is 

L 2 C A 



U(R) 



2R 2 R 



R 2 



(13) 



(14) 



We will see below that the three parameters (L, C, A) can 
be effectively reduced to two. 



B. Specifying the initial-value problem 



Solutions of (flU)) and (|12b[) are specified by initial con- 
ditions (R,R,ip,ip)(to) = (Ro,Vo,tpo,ujo) a t the initial 
time t . The discussion of the dynamical behavior of R 
is most effectively done in terms of the effective poten- 
tial. Moreover, since perturbations are best discussed in 
terms of dimensionless parameters, we also introduce a 
length scale and a time scale that appropriately charac- 
terize the dynamical perturbation and the solution to be 
perturbed. 

The length scale is defined as the radius at which the 
acceleration due to the cosmological expansion has the 
same magnitude as the two-body attraction. This hap- 
pens precisely at the critical radius 



Rr 



C_ 



1/3 



(15) 



For R < R c the two-body attraction dominates, whereas 
for R> R c the effect of the cosmological expansion is the 
dominant one. 

In order to gain an understanding of the length scales 
of the critical radius it is instructive to express it in terms 
of the physical parameters. In the case of gravitational 
interaction we have C/\A\ = GM/(\q \HQ ) and thus 



Rr.= 



( RsR 2 H 
V %o| 



1/3 



(16) 



Inserting the approximate value go = —1/2 of the present 
epoch, this reduces to the Schucking radius (|4|). 

In the electromagnetic case, e.g. for an electron-proton 
system, we have C/|A| = {e 2 / m) / {\q \H^) . Defining, in 
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analogy with (J5]) , the length scale 
2e 2 



5.64- lCT 15 ni, 



the critical radius (fTS")) becomes 



R, 



R e R 2 H 



1/3 



30 AU. 



(17) 



(18) 



where in the last step we inserted qo = —1/2. This is 
about as big as the Neptune orbit! 

From (fTtj)) and (fTS)) one sees that, in both cases, a larger 
(smaller) |go| implies a smaller (larger) critical radius, 
according to expectations. 

So much for the length scale. The time scale is defined 
to be the period of the unperturbed Kepler orbit (a solu- 
tion to the above problem for A = 0) of semi-major axis 
i?o- By Kepler's third law it is given by 



7\ := 2;r I 



(19) 



It is convenient to introduce two dimensionless param- 
eters which essentially encode the initial conditions Rq 
and ujq. 



A 



:- (■ 



\2tt/T k 



L 2 



CRq 



a := signtA) ( £ =^ 



(20) 
(21) 



For close to Keplerian orbits A is close to one. For reason- 
ably sized orbits a is close to zero. For example, in the 
Solar System, where R < 100 AU, one has |or| < 10 16 . 
For an atom whose radius is smaller than 10 Bohr-radii 
we have |a| < 10~ 57 . 
Now, defining 

x(t):=R(t)/R , (22) 
equations (fT3|) and (|12b|) can be written as 

+ (2n/T K ) 2 u X:a {x) =e (23) 

x 2 <p = loq , (24) 

where e := E/rfi now plays the role of the energy- 
constant and where the reduced two-parameter effective 
potential u\ >a is given by 



A 



1 



u\, a {x) := —2 -x . 

2ar x 2 

The initial conditions now read 

(x,x,<p,(p)(t ) = {l,V /Ro,ipo,oj ) 



(25) 



(26) 



The point of introducing the dimensionless variables is 
that the three initial parameters (L, C, ^4) of the effective 
potential could be reduced to two: A and a. This will be 
convenient in the discussion of the potential. 




FIG. 1 The figure shows the effective potential u\ tCt for cir- 
cular orbits (for which A = 1 — a) for some values of a. The 
initial conditions are x = 1 and x — (see (|22|l '). At x — 1 the 
potential has an extremum, which for a < 1/4 is a local mini- 
mum corresponding to stable circular orbits. For 1/4 < a < 1 
these become unstable. The value a — corresponds to the 
Newtonian case. 



C. Discussion of the reduced effective potential 

Circular orbits correspond to extrema of the effective 
potential (fl"4"|) . Expressed in terms of the dimensionless 
variables this is equivalent to u' x a (l) = — A + 1 — a = 0. 
By its very definition (|20[) , A is always non- negative, im- 
plying a < 1. For negative a (decelerating case) this is 
always satisfied. On the contrary, for positive a (acceler- 
ating case), this implies, in view of (|2"Tj) . the existence of 
a critical radius, given by R Cl beyond which no circular 
orbit exists. These orbits are stable if the considered ex- 
tremum is a true minimum, i.e. if the second derivative 
of the potential evaluated at the critical value is positive. 
Now, u'l (1) = 3A — 2 — a = 1 — 4a, showing stability for 
a < 1/4 and instability for a > 1/4. For the accelerating 
case, in view of (f2"Tj) , this implies that the circular orbits 
are stable iff Rq is smaller than the critical value 



Rs, 



:= (l/4) 1 / 3 i? B «0.63.R c 



(27) 



where 'sco' stands for 'stable circular orbits'. 

Summarizing, we have the following situation: in the 
decelerating case (i.e. for negative a or, equivalently, for 
negative ^4) stable circular orbits exist for every radius 
i?o; one just has to increase the angular velocity by some 
amount stated below in (f2"8"]) . On the contrary, in the 
accelerating case (i.e. for positive a, or, equivalently, for 
positive A), we have three regions: 

• i?o < -Rsco, where circular orbits exist and are sta- 
ble, 
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• Rsco < Rq < Rc-, where circular orbits exist but are 
unstable, and 

• Rq > R c , where no circular orbits exist. 

Generally, there exist no bounded orbits that extend be- 
yond the critical radius i? c , the reason being simply that 
there is no R > R c where U'(R) > 0. Bigger systems will 
just be slowly pulled apart by the cosmological acceler- 
ation and approximately move with the Hubble flow at 
later times. 9 Modifications of this strict qualitative dis- 
tinction impli ed by time dependencies o f A in (|14p were 
discussed in ( Faraoni and Jacques! . l2007h . 

Turning back to the case of circular orbits, we now 
express the condition for an extrema derived above, A = 
1 — a, in terms of the physical quantities, which leads to 



uj a = {2n/T K )^/l - sign(A)(i? /Rc 



(28) 



This equation says that, in order to get a circular or- 
bit, our planet, or electron, must have a smaller or big- 
ger angular velocity according to the universe expanding 
in an accelerating or decelerating fashion, respectively. 
This is just what one would expect, since the effect of 
a cosmological 'pulling apart' or 'pushing together' must 
be compensated by a smaller or larger centrifugal forces 
respectively, as compared to the Keplerian case. Equa- 
tion (|28| represents a modification of the third Kepler 
law due to the cosmological expansion. In principle this 
is measurable, but it is an effect of order (Rq/R c ) 3 and 
hence very small indeed; e.g. smaller than 10~ 17 for a 
planet in the Solar System. 

Instead of adjusting the initial angular velocity as 
in (|28|) . we can ask how one has to modify ro in or- 
der to get a circular orbit with the angular velocity 
ll>q = 2-k/Tk- This is equivalent to searching the min- 
imum of the effective potential (|25[) for A = 1. This con- 
dition leads to the fourth order equation ax A — x + 1 = 
with respect to x. Its solutions can be exactly writ- 
ten down using Ferrari's formula, though this is not il- 
luminating. For our purposes it is more convenient to 
solve it approximately, treating a as a small perturba- 
tion. Inserting the ansatz x m i n = cq + c\a + 0(a 2 ) 
we get Co = C\ = 1. This is really a minimum since 
u'{ a (x m { n ) = 1 + O(a) > 0. Hence we have 



J?min — Ro 



(29) 

This tells us that in the accelerating (decelerating) case 
the radii of the circular orbits with luq = 2t:/Tx be- 
comes bigger (smaller), again according to physical ex- 
pectation. As an example, the deviation in the radius 
for an hypothetical spacecraft orbiting around the Sun 



at 100 AU would be just of the order of 1 mm. Since it 
grows with the fourth power of the distance, the devia- 
tion at 1000 AU would be of the order of 10 meters. 



IV. GENERAL-RELATIVISTIC TREATMENT FOR 
ELECTROMAGNETICALLY-BOUNDED SYSTEMS 

In this section we show how to arrive at Q12[) from a 
relativistic treatment of an electromagnctically bounded 
two-body system embedded (without back-reaction) into 
an expanding (spatially flat) universe. This implies 
solving Maxwell's equations in the cosmological back- 
ground ([50)1 for an electric point charge (the proton) 
and then integrate the Lorentz equations for th e motion 
of a particle (electron) in a bound orbit (cf. ([Bonnorl . 
1999)). Equation p2[) then appears in an appropri- 
ate slow-motion limit. However, in order to relate this 
straightforward method to a famous argument of Dicke & 
Peebles, we shall proceed by taking a slight detour which 
makes use of the conformal properties of Maxwell's equa- 
tions. 



A. The argument of Dicke and Peebles 



In reference ( Dicke and Peebles! I1964T) Dicke & Pee- 
bles presented an apparently very general and elegant 
argument that purports to show the insignificance of any 
dynamical effect of cosmological expansion on a local sys- 
tem that is either bound by electromagnetic or gravita- 
tional forces and which should hold true at any scale. 
Their argument involves a rescaling of spacetime coor- 
dinates, (t, x) i ► (Xt, Ax) and certain assumptions on 
how other physical quantities, most prominently mass, 
behave under such scaling transformations. For exam- 
ple, they assume mass to transform like m i— > A _1 m. 
However, their argument is really independent of such 
assumptions, as we shall show below. We work from first 
principles to clearly display all assumptions made. 

We consider the motion of a charged point particle in 
an electromagnetic field. The whole system, i.e. particle 
plus electromagnetic field, is placed into a cosmological 
FLRW-spacetime with flat (k = 0) spatial geometry. The 
spacetime metric reads 



g = c 2 dt 2 -a 2 (t)(dr 2 



r 2 g s < 



)■ 



(30) 



where 



g s 2 = d9 2 + sin 2 6 dip 2 (31) 

denotes the metric on the unit two-sphere in standard 
coordinates. We introduce conformal time, t Cl via 

rt dt' 



to = fit) 



a(t>) 



(32) 



9 This genuine non-perturbative b ehavior was not seen in th e per- 
turbation analysis performed in (Coopcrsto ck et aU [l998). 



by means of which we can write (f30j) in a conformally flat 
form 

g = a 2 c (t c )(c 2 dt 2 c -dr 2 -r 2 g S 2)=a 2 c (t c )r 1 , (33) 
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where tj denotes the flat Minkowski metric. Here we 
wrote a c to indicate that we now expressed the expansion 
parameter a as function of t c rather than t, i.e. 



a c := ao f 1 



(34) 



The electromagnetic field is characterized by the tensor 
F^, comprising electric and magnetic fields: 



p _ ( E n /c 
^ v ~ \ -E lr —f R 



In terms of the electromagnetic four-vector potential, 
= (cp/c,—A), one has 

Fjxv = d„A v - d v A„ = V„A V - V V A„ , (36) 

so that, as usual, E = — V0 — A. The expression for the 
four-vector of the Lorentz-force of a particle of charge 
e moving in the field is eF^ L v u u , where w M is the 
particle's four velocity. 

The equations of motion for the system Particle + EM- 
Field follow from an action which is the sum of the action 
of the particle, the action for its interaction with the 
electromagnetic field, and the action for the free field, all 
placed in the background ([30)1 . Hence we write: 



S — Sp + Sj + Sf 



(37) 



where 
S P 
Si 



2 I dr = -mc I yj g{z' , z') d\ , (38a) 
L ,/.(•'' = - , / A Al (z(A))z'^ l dA 



d 4 xA fl (x) / dX e 8^ (x - z{X)) z'», 



S F = 



-1 



g <rg vf> F^Fcf 
d i x rf a if^ F^ v F a p ■ 



(38b) 



(38c) 



Here A is an arbitrary parameter along the worldline z : 
A i > z(X) of the particle, and z' the derivative dz/dX. 
The differential of the proper time along this worldline is 



dr 



y /gJ^r)dX = yJg^(z(X))^^d\. (39) 



It is now important to note that 1) the background met- 
ric g does not enter (|38b[) and that (|38c[) is conformally 
invariant (in 4 spacetime dimensions only!). Hence the 
expansion factor, a(t c ), does not enter these two expres- 
sions. For this reason we could write (|38c[) in terms of 
the flat Minkowski metric, though it should be kept in 
mind that the time coordinate is now given by conformal 
time t c . This is not the time read by standard clocks 
that move with the cosmological observers, which rather 



show the cosmological time t (which is the proper time 
along the geodesic flow of the observer field d/dt). 

The situation is rather different for the action (|38a|) 
of the particle. Its variational derivative with respect to 
z(X) is 



8S V 



Szt*(X) 



d 

dX 



(35) 

We now introduce the conformal proper time, r c , via 



(40) 



dr c = (1/c) y/r)(z',z')d\ = (1/ca) y/g(z' , z>) dX . (41) 

We denote differentiati on with respect to t c by an over- 
dot, so that e.g. z'/ y/g(z', z') = z/ca. Using this to 
replace z' by zy/ ' g(z' , z')/ca and also g by a 2r q in (|40|) 
gives 



5S„ 



6z»{X) 
where we set 



ac 



{^ a z a -P^4>,a} (42) 



a=:exp(0/c 2 ) and P« := ^ 



(43) 



Recalling that SS P = J j^Sz^dX = J j^Sz^dr. 
and using (|4"Tj) , (|4^|) is equivalent to 



5S„ 



o~z»{t c ) 



ma(z a - P«<t>, a ) , 



(44) 



where from now on we agree to raise and lower indices us- 
ing the Minkowski metric, i.e. rj^ v = diag(l, — 1, — 1, — 1) 
in Minkowski inertial coordinates. 

Writing (|38b[) in terms of the conformal proper time 
and taking the variational derivative with respect to z(t c ) 
leads to 8Si / 8z^{t c ) = —eF^ a z a , so that 



SS 



5z»(t c ) 



i(z^P^^ a )-eF^ a z a . (45) 



The variational derivative of the action with respect to 
the vector potential A is 



SS 



= d a F" a (x)-e dT c 6{x-z(T c ))z»(T c ). (46) 



SA M (x) 



Equations (|45j) and (|46|) show that the fully dynamical 
problem can be treated as if it were situated in static flat 
space. The field equations that follow from (j46|) are just 
the same as in Minkowski space. Hence we can calcu- 
late the Coulomb field as usual. On the other hand, the 
equations of motion receive two changes from the cos- 
mological expansion term: the first is that the mass m 
is now multiplied with the (time-dependent!) scale fac- 
tor a, the second is an additional scalar force induced by 
a. Note that all spacetime dependent functions on the 
right hand side are to be evaluated at the particle's loca- 
tion z(t c ), whose fourth component corresponds to ct c . 
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Hence, writing out all arguments and taking into account 
that the time coordinate is t c , we have for the equation 
of motion 

F»(z)z a 



physical coordinates. 10 The only additional term next 
to y that survives is precisely the familiar acceleration 
term ([7]). Inserting (|4"9")l into l|48p. whose right-hand 
side we now specialize to a pure electric Coulomb field, 
E(z) — Qz/\z\ 3 and B(z) = 0, we arrive at 



u. —♦/••/ \ eQ _ 
y-y (a/a) = — —y. 

m \y\ 



After introducing polar coordinates in the orbital plane 
we exactly get (|T2|) . 



ma c (z° I c) 

- (-cV" + z"z a )d a \na c {z°/c) (47a) 

= ma c {z°/c) F ^ za 

- (-cY° + z^z°/c) a' c (z°/c)/a c {z Q /c) , (47b) 

where a' c is the derivative of a c . 

So far no approximations were made. Now we write 
z^ = 7(c,u), where v is the derivative of z with respect B Exact con dition for non-expanding circular orbits 

to the conformal time t c , henceforth denoted by a prime, 
and 7 = 1/yl — v 2 /c 2 . Then we specialize to slow mo- 
tions, i.e. neglect effects of quadratic or higher powers in 
v/c (special relativistic effects). For the spatial part of 
(|47b|) we get 



(50) 



+ z' (a'Ja c ) = (E + z'xB), 



(48) 



where we once more recall that the spatial coordinates 
used here are the comoving (i.e. conformal) ones and the 
electric and magnetic fields are evaluated at the particle's 
position z(t c ). 

From the above equation we see that the effect of cos- 
mological expansion in the conformal coordinates shows 
up in two ways: first in a time dependence of the mass 
which scales with a c , and, second, in the presence a fric- 
tion term. Dicke & Peebles neglect the friction term 
and simply conclude as follows: In the adiabatic approx- 
imation, which is justified if typical time scales of the 
problem at hand are short compared to the world-age 
(corresponding to small e 2 in (|98b[l ). the time-dependent 
mass term leads to a time varying radius in comoving 
(or conformal) coordinates of r{t c ) cx l/a c (t c ). Hence 
the physical radius (given by the cosmological geodesi- 
cally spatial distance), r* = a c r, stays constant in this 
approximation. Hence, within this approximation, elec- 
tromagnetically bound systems do not feel any effect of 
cosmological expansion. 

But what does 'this approximation' refer to? We will 
see that it amounts to neglecting precisely the leading 
order contributions. This is easy to see if we cast (|48p 
into physical coordinates, given by the cosmological time 
t and the cosmological geodesic spatial distance r* := 
a(t)r. We have dt c /dt = 1/a and the spatial geodesic 
coordinates are y := a(t)z. Denoting by an overdot the 
time derivative with respect to t, the two terms on the 
left hand side of (flBl become 



z"(a'Ja c ) 



a(yH 2 - $H) + aQ - yd/a) , (49a) 
-a(yH 2 - jpT) , (49b) 



where H = a/a. This shows that the friction term can- 
cels against the first-order derivative terms in y and a 
that one gets in rewriting the left-hand side of (|48p in 



In dBonnorl . [l999) a necessary and sufficient condi- 
tion for the existence of non-expanding orbits is derived 
for the electron-proton system in a spatially flat FLRW 
spacetime. Here 'non-expanding' is defined as of con- 
stant areal radius. This condition follows directly from 
the Lorentz equation of motion for the electron in the 
external electric field of the proton, the normalization 
condition of the electron's four-velocity, and the condi- 
tion of constancy of the areal radius. In our notation, in- 
troducing the dimcnsionless quantities h(t) := RH(t)/c, 
I := L/Rc, and /Lt := R e /2R, the conditions for the exis- 
tence of non-expanding circular orbits reads as follows: 



R. (1-h 2 ) 3 / 2 ( l 2 + h 2 

7 - (i+/2)i/2 r~ 



(51) 



Recall that R e is defined in (fT7|) and H(t) and L denote 
the Hubble function and, respectively, the (conserved) 
electron's angular momentum per unit mass. The above 
condition is a first-order autonomous ODE for the func- 
tion h(t), and hence for the Hubble function H(t). This 
is the constraint on the spacetime (more precisely, on 
the scale factor a(i)) that one gets by imposing the ex- 
istence of non-expanding circular orbits for two oppo- 
sitely charged point masses. If such orbits exist, (|5T|) 
amounts to the generalization of Kepler's third law to 
FLRW spacetimes, which here gives a relation between 
the scale function on one hand and the orbital parame- 
ters R and L as well as the field-strength parameter R e 
on the other. Recall that in Newtonian physics the third 
Kepler law is, in our notation, simply given by I 2 = ji. 

The easiest solutions of (|5Tj) are of course the station- 
ary ones, that is with h(t) = ho, for some constant ho. 
This means that the scale factors is exponentially ex- 
panding, 



a(t) = a cxp(H Q t) , 



(52) 



10 Since the friction term cancels, the critical remark [27] in 
l|Adkins et all EqOT\ ) regarding its magnitude is based on a mis- 
understanding. 
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where Hq :— h$c/ R and ao is some positive constant. In 
other words, the spacetime is given by the de Sitter so- 
lution (A-dominated universe). In this case ([5"Tj) reduces 
to 



l 2 + h 2 



(53) 



Notice that a larger Hubble parameter, hence a larger 
ho, makes the l.h.s. larger. Consequently, (|5"3")) tells us 
that with a larger Hubble parameter we must give to the 
electron a smaller angular velocity (smaller I) in order to 
keep it on a non-expanding circular orbit with the same 
radius. This, according to intuition, is in order to com- 
pensate the extra cosmological pull with a reduced cen- 
trifugal term. In case of Minkowski spacetime (ho = 0) 
the above relation reads i 2 /vl + I 2 = /x, hence one can 
interpret the factor 1 / vl + P as a special- relativistic cor- 
rection to the Newtonian relation I 2 — /i. The largest ra- 
dius at which, in an FLRW spacetime with exponentially- 
growing scale factor, there is a non-expanding orbit fol- 
lows from ([53]) in the limit I — > 0. In this limit the condi- 
tion reduces to h^/y/l — h\ — [i, which, for small param- 
eters ho, simplifies to h 2 , sa /i. Solving for R this gives 
the radius (i? e i?| f /2) 1 / 3 , which, taking into account that 
qo = — 1 because of (j5"2")) . exactly corresponds to the crit- 
ical radius (fT8| . 

The other (non-stationary) solutions of (151|) can also be 
found. After separation of variables and an elementary 
integration one gets t as function of h in terms of trigono- 
metric functions composed with inverse hyperbolic func- 
tions. This exact expression is again not very illuminat- 
ing and cannot generally be explicitly inverted so as to 
obtain h in terms of elementary functions of t. How- 
ever, if we make use of the smallness of the parameters 
/i, I 2 , and h 2 , a leading-order expansion in these quanti- 
ties gives a much simpler expression for t(h) which can 
be explicitly inverted. In fact, this approximate solution 
h(t) is obtained much quicker by solving ([51]) with the 
right-hand side being replaced with its leading-order ex- 
pansion in the mentioned quantities, that is, by solving 



R; 



V-l 2 



(54) 



Here /j — I 2 is a constant which depends on the orbit pa- 
rameters. One must now distinguish between three cases: 
(a) ft — I 2 =: k 2 > for some positive k, (b) /i — I 2 =: 
—v 2 < for some positive v, and (c) fi — I 2 = 0. Re- 
calling the Newtonian relation I 2 — /i, orbits in the three 
cases have an angular momentum which is, respectively, 
smaller, bigger, and equal to the Newtonian one. In- 
tegrating (|54"1) we get, putting w.l.o.g. t = 0, h(t) = 
k tanh(Kct/i?), h(t) — —vt&n{vct/R), and h(t) = R/ct, 
for the cases (a), (6), and (c), respectively. Then, inte- 
grating once and exponentiating the result, we get the 
corresponding scale functions: 

(a) Case /i — I 2 =: k 2 > (non-expanding orbits have 



sub- Newtonian angular momentum) 



net \ 



a(t) = a cosh ( — j , t G [0, oo) . (55a) 



(b) Case /i — I 2 =: —v 2 < (non-expanding orbits have 
super-Newtonian angular momentum) 



a(t) = a cos 



vet 
~R 



t € 



ttR 
2vc 



(55b) 



(c) Case fi — I 2 = (non-expanding orbits have Newto- 
nian angular momentum) 



a(t) = aot , 



t e (0,oo). 



(55c) 



In all three cases (a), (b), and (c) ao is a positive con- 
stant and the acceleration term a/a is a constant which is 
positive, negative, and zero, respectively. Hence, as one 
would intuitively expect, the non-expanding orbits have 
an angular momentum which is smaller, larger, or equal 
the Newtonian one, depending on whether the accelera- 
tion factor d/a is positive, negative, or zero. 



V. GENERAL-RELATIVISTIC TREATMENT FOR 
GRAVITATIONALLY-BOUNDED SYSTEMS 



As advertised in Section III.B1 we now wish to dis- 
cuss exact solutions that may represent quasi-isolated 
spherically-symmetric gravitating systems 'embedded' 
into cosmological spacetimes. As regards the meaning 
of 'embedded' we distinguish between the strategies of 
'matching' and 'melting', as outlined in Sect. III. Bl 



A. Spherically-symmetric matchings 

The complexity and non-linearity of Einstein's equa- 
tions make it a very difficult task to construct a suitable 
variety of exact solutions which serve as realistic models 
for actual physical situations. Often exact solutions are 
only known for highly idealized situations, typically with 
high degrees of symmetry, in which the field equations 
sufficiently simplify. One way to construct new solutions 
(in a suitable sense, see below) from old ones is to glue 
them across suitably chosen hypersurfaces along which 
the matter distribution may become singular due to sur- 
face layers. This approach was pioneered by Lanczos in 
the early 1920s ( Lanczos|. 1924 ) and p ut into geom etric 



form by Darmo is (iDarmois 
also S 21.13 of (Misner et al 



19271) a nd (|lsrae]Lll96fj) : see 
19731 ). In this section, 



der the assumption of spherically symmetry, we present 
a new alternative set of conditions which are equivalent 
to the old ones. The new conditions only involve scalar 
quantities, are easy to verify, and have good ph ysical in- 
terpr etations. More details are contained in ( Carreral 
I2009D . 
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Here we shall restrict to piecewise continuous matter 
distributions without singular (^-distribution like) sur- 
face layers, as e.g. in the presence of stars with sharply 
defined surfaces. Einstein's equations can the be satisfied 
for piecewise twice continuously differentiable fields, if 
the field equations at the location of the matching hyper- 
surface are replaced by their one-dimensional e-interval 
integrals in normal direction to the hypersurface. The 
condition that two twice continuously differentiable so- 
lutions (in the ordinary sense) can be matched into a 
piecewise twice continuously differentiable solution (in 
the re-interpreted sense just explained) is then simply 
given by the so-called 

Darmois junction conditions (DJC). For a non- 
null matching hypersurface r , (i) the induced metric gr 
and (ii) the extrinsic curvature Kr shall be continuous 
through r . 

Let us pause for a moment to say a few more words 
about the notion of 'continuity through r\ Gluing to- 
gether two pieces of spacetimes means the following: 
Initially one has two spacetimes, say (Ai + ,g + ) and 
(M~ ,g~), with oriented boundaries r + and r~ , respec- 
tively. Given a diffeomorphism (f> : <j> : r + — > r~~ between 
the boundaries, the glued spacetime is the quotient of the 
disjoint union of M + and M.~ under the identification 
of each point of p G r + with <f>(p) S r~. The match- 
ing hypersurface r is now the common image of r + and 
r~ after identification in the quotient spacetime. Now, 
a tensor field T is said to be continuous through r if 
T\p+ equals T\p- under the push-forward action of the 
diffeomorphism (f>, hence if </>*(T| r +) = T\ r -. 

Let us now return to the DJC and, in particular, their 
physical interpretation. If n is a continuous choice of 
unit normal of _T, it implies that 



T(n, •) is continuous through r . 



(56) 



This follows directly from the expressions (|C8[) of the 
Einstein tensor given in Appendix [C] If r is time- 
like and hence n spacelike, (I5"6")) just states the continu- 
ity of the normal components for the energy-momentum 
flux-densities, whereas their tangential components to- 
gether with the energy density may jump across r. In 
the absence of surface layers this continuity condition 
is just a physically obvious consequence of local energy- 
momentum conservation, whereas jumps in, say, the 
energy-density must clearly be allowed for. For com- 
pleteness we note that for a spacelike matching surface 
(|S"6")) states the continuity of the densities of energy and 
momentum as measured by an observer moving along n 
(taken to be future pointing) , whereas the corresponding 
currents may jump. 

Let now restrict our attention to spherical symmetric 
spacetimes glued along hypersurfaces of spherical symme- 
try. This means that the latter are left invariant, as set, 
under the action of SO (3). We recall that the structure 
of a spherically symmetric spacetime is that of a warped 
product M. — B x R S 2 of a two-dimensional Lorentzian 



manifold B and the two-sphere by means of the warp- 
ing function, R : B — ► R+, called the areal radius. The 
matching hypersurfaces are then of the form r = -fXuS 2 , 
where 7 is the projection of r into B and is called the 
matching curve. The DJC should then reduce to appro- 
priate conditions along the curve 7. Indeed, Theorem [T] 
below shows that, in the spherically symmetric case, the 
DJC are equivalent to the following 

Spherically symmetric junction conditions 
(SSJC). Let r be a smooth, non-null, spherically sym- 
metric matching hypersurface between two spherically 
symmetric spacetimes and n a continuous choice of unit 
normal vector field on T '. Denote with 7 the projection 
of r onto B. Moreover, let v the (unique up to a sign) 
spherically symmetric, unit vector field on r orthogonal 
to n. The following four functions 
(i) the arc-length 0/7, 

(ii) the extrinsic curvature 0/7 in B: g(n, V„i>), 
(Hi) the areal radius R, 
(iv) the Misner-Sharp energy E, 
shall be continuous through the matching curve 7. 



Th e Mis ner-Sharp (MS) energy ( Misner and Sharp! . 
Il964h and ( Hernandez and Misner . [l966h . which we will 
denote by E, is a concept of quasi-local mass that can 
be defined in presence of spherical symmetry, and which 
proves useful for computational and interpretational pur- 
poses. It is a function defined in purely geometrical man- 
ner as follows. Given a point p of spacetime, compute 
the sectional curvature K of the plane tangent to the 
two-dimensional S f O(3)-orbit through p and multiply this 
with minus 11 one-half of the third power of the areal ra- 
dius: 



E := -±R 3 K. 
From (|D9|) we immediately read off 

E = ^(l + {dR,dR)) 



(57) 



(58) 



which provides a convenient expression for the computa- 
tion of the MS energy. We shall show in Appendix iDl that 
the MS energy is the charge of a conserved current and 
how it depends on the energy-momentum tensor for the 
matter. There we will also briefly discuss its Newtonian 
limit. This allows to interpret it as amount of active grav- 
itational energy contained in the interior of the sphere of 
symmetry (S f O(3)-orbit) through p. There we also show 
that the MS energy at p is equal to the Hawking quasi- 
local mass of the two-sphere of symmetry (S*0(3)-orbit) 
through p and hence converges to the Bondi-mass at null 
infinity and, in an asymptotically flat spacetime, to the 
ADM mass at spatial infinity (for the latter two issues 



The minus sign here is just a relict of our signature choice. 
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see (|Szabadost 12001 and also (|Havwardl Il996t )). More- 
over, we will give the decomposition of the MS energy in 
its Ricci and Weyl parts; see (|D34|) . 

The name 'Misner-Sharp energy' seems now to be 
established in the literature, however one should sa~ 



that this mass concept goes back at least to (|Lemaitr< 
I1933TI 12 . which gives a coordinate expression for 
it. Its geometric definition (|5T[) was first given 
in ([Hernandez and Misneii Il966l ) and its interpretation 
as the charge of a conserved current was first derived 
in (lKodama|. [19 80). Later, an alternative definition was 
Tl990f ): There it is showed that the MS 



given m 

energy can defined in terms of the norm of the Killing 
fields generating the isometry group SO(3), leading di- 
rectly to (l58l) . Fu rther releva n t studi es of the MS en- 
ergy ar e dCahill and McVittieL 1970a l a) an d, more re- 



cently, (|Burnettl . 1 19911: iHavwardl ll996[Tl998l ). 

Some comment are needed on the above SSJC. First we 
note that, since n and v are spherically symmetric and 
hence tangent to B we have in view of (|D3|) and (|C6|) 
that g(n, V v v) = g B {n B , B V VB v B ) = e(n B ) K*(v B , v B ). 
Hence, the quantity in (ii) is indeed (up to a possible 
sign) the extrinsic curvature of the curve 7 in B. Second, 
note that this quantity, being quadratic in v, does not de- 
pend on the sign choice of v. Third, since the matching 
hypersurface has the structure 7 x S* 2 , the words 'contin- 
uous through the curve 7' can be interchanged with the 
words 'continuous through the hypersurface f", depend- 
ing on ones preference to think four- or two-dimensional. 

A great advantage of the SSJC is that they are very 
easy to verify: one simply has to impose continuity 
on four scalar s along the matching curve in the two- 
dimensional base manifold B. Dealing with scalars, since 
their value is independent on the particular coordinate 
choice, one does not need to worry about introducing 
new coordinates in both spacetimes to be glued, in order 
to get the different metrics in a form which is comparable. 
This is indeed an ingrate task: in general, these coordi- 
nates are only needed in order to check if the junction 
conditions are satisfied, and for nothing more. In pres- 
ence of spherical symmetry all this can be circumvented 
by using our new junction conditions. 

Furthermore, the SSJC have a good physical interpre- 
tation: The continuity condition of both, the areal radius 
as well as of the MS energy, can be read as equilibrium 
condition for the gravitational pull acting from opposite 
directions onto (fictitious) test masses at the location of 
the matching surface. Concerning the continuity of the 
extrinsic curvature of the matching curve, we note the 
following: In the case where the matching hypersurface 
r is timelike, let 7 = n(T) be a timelike curve in B and 
v is future-pointing tangent. One can think of 7 as the 
'matching observer's' worldline. Hence, in the timelike 
case, the extrinsic curvature is nothing but the acceler- 



ation of the matching observer. In the spatial case, on 
the contrary, n is timclike, 7 is spacelike (and v tangent 
to it). One can choose n to be future-pointing and then 
think of it as an observer field defined along the spatial 
(1-dimensional) slice 7. In view of (|C7[) one then sees 
that the extrinsic curvature of 7 is exactly 13 the shear- 
expansion of n in direction of v 'radial direction'. 
In Appendix [B] we prove the following 

Theorem 1 (Equivalence of the junction conditions). 
Let r be a smooth, non-null, spherically symmetric 
matching hypersurface between two spherically symmet- 
ric spacetimes and n a continuous choice of unit normal 
vector field on F. Assume, moreover, that the areal radii 
of the two spacetimes are C 1 functions in an open neigh- 
borhood of the matching hyper surf aces. Then the DJC 
are equivalent to the SSJC. 

Now let us suppose we are faced with the following 
situation: We are given two spherically symmetric solu- 
tions of Einstein equation and we want to know if they 
can be matched together at all, and if so, how to char- 
acterize the curve 7 (respectively, the hypersurface T) 
along which this is possible. Answers to these questions 
will be provided by the junction conditions SSJC. Note 
that a timelike or spacelike curve in the two-dimensional 
base manifold (B,g B ) can be described simply by a func- 
tion R(t), where R is the areal radius and r the curve's 
arc-length which, in the timelike case, corresponds to the 
matching observer's proper time. The conditions (i) and 
(Hi) of the SSJC are then equivalent to the condition 
that the functional dependence R(t) must be the same 
(up to a trivial translation in r) in both spacetimes to be 
matched. 



B. The Eisenstaedt theorem 

Perhaps the simplest attempt to model a compact 
body (star) in an expanding universe is trying to inglo- 
bate it in a FLRW spacetime and to assume, for sim- 
plicity, that the body is spherically symmetric. A direct 
consequence of the SSJC is the fo llowing intui t ive ap peal- 
ing theorem due to Eisenstaedt (jEisenstaedtl . I1977T ): 

Theorem 2 (Eisenstaedt, 1977). Excise the full world- 
tube W ro of a comoving ball of comoving radius r^ from a 
FLRW spacetime and insert instead a spherical symmet- 
ric inhomogeneity (hence a piece of a spherically sym- 
metric spacetime together with a related matter model, 
satisfying Einstein's equations). Then a necessary condi- 
tion for the resulting spacetime to satisfy Einstein's equa- 
tion is that the MS energy of the inserted inhomogeneity 
equals that of the excised ball. 



12 For an English translation see llLemaitrd Il997t) . 



13 Recall that, because of our signature choice, the restriction of 
the metric to spacelike directions is negative definite. 
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This says that the mean energy density (measured with 
the MS energy) of spherically-symmetric inhomogcncitics 
must be the same as the one of the FLRW spacetime. 
That the Eisenstaedt Theorem is a consequence of the 
Darmois junct ion conditions was already pointed out 
in (|Hartl 120061 ). 



C. The Einstein-Straus vacuole revisited 

As another application of the above described match- 
ing procedure we revis i t the Ein s tein-Strau s solu - 
tion jEinstein and Strausl Il945l Il946h ISchuckinel . Il954 . 
which originally consists on a Schwarzschild spacetime 
(called 'vacuole') matched to a dust FLRW universe with 
zero cosmological constant. Later, this model was gen- 
eralized also to the case of a non- vanishing cosmologi- 
cal constant ( Balbinot et all Il988h . We treat here the 



general case of an arbitrary cosmological constant and 
show that the SSJC allow substantial simplifications of 
the computations. This technique can also be applied to 
Bonnor's vacuole construction in LTB spacetimes. 

Notice that the matching condition (JS5J) implies, in 
particular, that the pressure must be continuous through 
the matching hypersurface. Since the interior is a vacuum 
spacetime, it follows that the pressure must vanishes also 
on the exterior part of the matching hypersurface and 
hence everywhere on the FLRW spacetime. That is why 
one has to restrict to dust FLRW spacetimes. 

Since we leave the cosmological constant A arbitrary 
(it may be positive, negative, or zero) the inner re- 
gion is given, respectively, by a Schwarzschild-de Sitter, 
Schwarzschild-anti-de Sitter 14 , or Schwarzschild space- 
time (all abbreviated henceforth by SdS). Recall that the 
SdS spacetime is given by the vacuum solution to Ein- 
stein equation with cosmological term 



g SdS = V(R) dT 2 - V{R)- l dR 2 - R 2 g S 2 



(59a) 



together with the matter energy-momentum tensor 15 
T = gu®u, where u = d/dt is the (geodesic) velocity 
field of the cosmological dust and g is the matter energy 
density, which depends on t only. Here r is the comov- 
ing radial coordinate and k is a constant which takes the 
values 0, —1, +1, depending on whether the spatial slices 
have zero, negative, or positive curvature, respectively. 
The Einstein equation is then equivalent to the following 
set given by the Friedmann equation and a 'conservation 
equation': 



ga — constant, 



(61a) 
(61b) 



where the constant C := 8ng ag/3 depends on the initial 
conditions a := a(t ) and g„ :— g(t ) at some 'initial' 
time t . Here, the dot denotes differentiation with respect 
to t or, which is the same, along u. 

The central question is now the following: How shall we 



cut hypersurfaces r s , 



7sd 



x S 2 and r plRW — 7flrw x 



S in the spacetimes SdS and, respectively, FLRW in 
order that the resulting pieces can be matched? In order 
to apply the SSJC we have to compute the MS energy 
for both spacetimes. For the FLRW spacetime one has 
-Rflrw(£, r) — a(t)r and hence dR FllIfw = ardt + adr and 
(dR, e£i?) FLRW = g^R^R^ = r 2 (k + a 2 ) - 1. From the 
definition (j58"| one then gets, using Friedmann's equation, 



Air 3 . 

y fi FLEw(? + QA) 



(62) 



where g\ :— A/87T is the energy density associated with 
the cosmological constant. Notice that this expression, 
as well its derivation, is completely independent of the 
specific equation of state of the fluid and does not depend 
on the spatial curvature k. For the SdS case one has 
R SdS = R and thus (dR, dR) SdS = — V(R) and 



where 



47F n3 



(63) 



(59b) 



Above, g s -2 denotes the metric on the unit two- 
sphere (f3"Tj) and to is a constant which represents the 
central mass. 

A dust FLRW spacetime is given by the metric 



Now, the last two conditions of the SSJC, that is 
the continuity of areal radius and MS energy across 
the matching hypersurface (yet to be determined), are 
equivalent to the continuity of the areal radius -R F lrw = 
R SdS —: R, together with the suggestive relation 



47T 



-R 3 g. 



(64) 



This two conditions already determine the matching hy- 
persurface. Indeed, inserting R = R FllIiV j = a(t)r in 



The Schwarzschild-(anti-) de Sitter metric (|59afl is often called 
the Kottler solution, after F riedrich Kottle r, who was the first 
to write down this metric in ijKottlcr. 1918). M ore details on its 
analytic and global structure may be found in l lGeverl . Il98dt) . 



Throughout we denote the metric-dual (1-form) of a vector u 
by underlining it, that is, u:= g(u, ■ ) is the 1-form metric-dual 
to vector u. In local coordinates we have u = rt M 9 M and u = 
Uf+dx^, where := g l _ lv u v . 



17 



and using the relation g(t)a 3 (t) — const, valid for dust 
FLRW models, one obtains the matching radius in terms 
of the FLRW comoving radial coordinate: 



ra := 



1/3 



constant . 



(65) 



Here a := a(t ), and similarly for p, where t is some 
fixed 'initial' time. This means that the matching ob- 
server moves, in the FLRW spacetime, along the inte- 
gral curve of it = d/dt with initial condition (t ,r ) and 
hence is comoving with the cosmological matter. 

So far we used the last two of the SSJC. As discussed 
above, the continuity of the areal radius and the arc- 
length (the proper time, in the timelike case) of the 
matching curve are equivalent to the equality of the func- 
tional dependencies R{t) (up to a possible trivial trans- 
lation in r) which describe the matching curves in the 
two spacetimes to be matched. Now, because of (|65|). 
the matching curve (worldlinc) in the FLRW spacetime 
is simply 



R(t) = a(r)r , 



(66) 



where a is the (unique) solution of the Friedmann equa- 
tion (|61a[) with initial condition a at t = t Q . (Recall 
that in FLRW the proper time of an observer moving 
along an integral line of u equals the cosmological time, 
hence r = t.) From what we said above, the same 
functional relation R(t) must hold also in the SdS — 
provided we identify R with the areal radius and r with 
the matching observer's proper time, both referred to the 
SdS spacetime. This determines the matching curve in 
SdS. 

Finally we need to show that the junction condition (n) 
is satisfied, hence that the matching observer's accelera- 
tions coincide. Looking at the matching worldline from 
the FLRW spacetime, it is immediately clear that it is 
geodesic, hence its acceleration vanishes. To conclude 
the matching procedure, we just have to check that this 
is also true for the matching worldline in the SdS space- 
time. For this, one has just to check that the function 
defined in (|66p satisfies the geodesic equation for a radial 
motion. The latter is given by 



R 2 + V(R) = 



(67) 



where e := g sds (d/dT, v) = constant and v = f 8/8T + 
R d/dR is the matching observer in the SdS space- 
time. (Equation (f6"T)) can be quickly derived from the fact 
that v(g sds (d/dT,v)) = 0, since d/dT is Killing and 
v geodesic. Inserting e := g SdS (d/dT,v) = V{R)f in 
the normalization condition 1 = g SdS (v,v) = V(R)T 2 — 
R 2 /V(R) one arrives immediately at (167p .) Now, insert- 
ing (|66p with (|65p in (|67p and using the Friedmann equa- 
tion (|61ap . one gets R 2 + V(R) = 1 — kr 2 . Hence, the 
geodesic equation (|S7)l is satisfied (with e 2 — 1 — kr 2 ) 
and herewith all the four junction conditions. 



D. The McVittie model 

Among all models discussed in the literature which rep- 
resent a quasi-isolated spherically-symmetric gravitating 
system melted into a cosmological spacetime, the one 
that is presumably best understood as regards its an- 
alytical structure as well as its p hysical assumptions is 
that of McVittie (iMcVittid. Il933h. tha nks to the careful 
analysis of Nolan (|NolanL I1998L Il999a0 bh. Here we shall 
restrict to the 'flat' or k = model, which interpolates 
between an exterior Schwarzschild solution, describing a 
local mass, and a spatially flat (i.e. k = 0) ambient FLRW 
universe. For simplicity we shall from now on refer to this 
model simply as the McVittie model. The cosmological 
constant is assumed to be zero, although this assumption 
is not essential (see the last paragraph of Section |VDT|. 

This is not to say that this model is to be taken at 
face value in all its aspects. Its problems lie in the re- 
gion very close to the central object, where the basic 
assumptions on the behavior of matter definitely turn 
unphysical. However, at radii much larger than (in geo- 
metric units) the central mass (to be defined below) the 
k = McVittie solution seems to provide a viable ap- 
proximation for the transition between a homogeneous 
cosmological spacetime and a localized mass immersed 
in it. We will now briefly discuss this model and look at 
its geodesic equations, showing that they reduces to (TT21) 
in an appropriate weak-field and slow-motion limit. This 
provides another and more solid justification for the New- 
tonian approach we carried out in Section [IIII 

The characterization of the McVittie model is made 
through two sets of a priori specifications. The first set 
concerns the metric (left side of Einstein's equations) and 
the second set the matter (right side of Einstein's equa- 
tions). The former consists in an ansatz for the metric, 
which can formally be described as follows: Write down 
the Schwarzschild metric for the mass parameter m in 
isotropic coordinates, add a conformal factor a 2 (t) to the 
spatial part, and allow the mass parameter m to depend 
on time. Hence the metric reads 



g = 



1 - m(t)/2r 



l + m(f)/2r 

4 



dt 2 



1 



in 



(*) 



2r 



(68) 



a 2 (t) (dr 2 +r 2 g 3 2). 



where g s 2 is given by pip . The metric (|68p is obviously 
spherically symmetric with the spheres of constant ra- 
dius r being the orbits of the rotation group. We will 
discuss below what this ansatz actually entails. For later 
convenience we also introduce the orthonormal tetrad 
{e^j-^o,--- ,3 with respect to (|68l) . where 16 



We write 



e M := Wd/dx^W^d/da 



V\9{v, v)\. 



(69) 
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and {x^} — {t, r, 9, ip}. 

The second set of specifications, concerning the matter, 
is as follows: The matter is a perfect fluid with density 
g and isotropic pressure p. Hence its energy-momentum 
tensor is given by 



T = QU®u+p(u^u — g) . 



(70) 



Furthermore, and this is where the two sets of specifi- 
cations make contact, the motion of the matter is given 
by 



u = e . 



(71) 



No further assumptions are made. In particular, an equa- 
tion of state, like p = p(g), is not assumed. The reason 
for this will become clear soon. 

Note that the vector field (|7Tj) is not geodesic for the 
metric (|6"5)) (unlike for the FLRW and Gautreau metrics), 
which immediately implies that the pressure cannot be 
constant. Being spherically symmetric, u is automati- 
cally vorticity free. The last property is manifest from 
its hypersurface orthogonality, which is immediate from 
(f6"8")l . Moreover, u is also shear free. This, too, can be 
immediately read off (|68|) once one takes into account 
the result that for spherically symmetric metrics van- 
ishing shear for a spherically symmetric vector field is 
equivalent to the corresponding spatial metric being con- 
formally related to a spherically symmetric flat metric. 
This is obviously the case here. 

The non-vanishing components of the Einstein tensor 
with respect to the orthonormal basis (|69p are: 

Ein(e ,e ) = 3F 2 , (72a) 

Ein( ei , ej ) = - (3F 2 + 2|f) <% , (72b) 

Ein(e ,ei) = J, (|) 2 (am)', (72c) 

where an overdot denotes differentiation along d/dt. Be- 
fore explaining the functions A, B, R, and F, we make 
the important observation that the Einstein tensor is spa- 
tially isotropic, where 'spatially' refers to the directions 
orthogonal to eo. By this we mean that Ein(ei, ej) oc 5ij 
or, expressed more geometrically, that the spatial restric- 
tion of the Einstein tensor is proportional to the spatial 
restriction of the metric. 

In (172")) and in the following we set: 

A(t,r) := l + m(t)/2r, B(t,r) := 1 - m(t)/2r , (73) 



and 



R(t,r) 



m(t) 
2r 



a(t) r , 



(74) 



where R is the areal radius for the McVittie ansatz 
and also 



F 



a 1 (am)' 

- + —f;- - 

a rB a 



(75) 



In passing we note that F has the geometric interpreta- 
tion of being one third the expansion of the vector field 
eo, that is, F — div(eo)/3. Hence (|72aj) could also be 
written in the form Ein(eo,eo) = (div(eo)) 2 /3. We will 
see later that the product am which appears in (|72c[) 
also has a geometric meaning: it is just the Weyl part of 
the MS energy; see (|8"Tj). 

Now, the non- vanishing components of the energy mo- 
mentum tensor (1701) with (1711) are: 



(76) 



T(e 0l e ) = g 1 T(e il e j )=pS i 



The (eo,ei) component of Einstein's equation therefore 
implies (am) = 0, which means that the Weyl part of the 
MS energy is constant. Physically this can be interpreted 
as saying that the central object does not accrete any 
energy from the ambient matter. Using the constancy of 
a m in (I75|) we immediately get: 



F=-=:H. 

a 



(77) 



Hence Einstein's equation is equivalent to the fol- 
lowing three relations between the four functions 
m(t), a(t), g(t, r), and p(t, r): 



(am) 
8irg 

8irp 



0. 









) 








a 


-»( 






a 



1 + m/2r 
1 - m/2r 



(78a) 
(78b) 

(78c) 



Note that here Einstein's equation has only three inde- 
pendent components (as opposed to four for a general 
spherically symmetric metric), which is a consequence of 
the fact, already stresses above, that the Einstein tensor 
for the McVittie ansatz (|68|) is spatially isotropic. 
Equation (|78ap can be immediately integrated: 



a(t) 



(79) 



where mo is an integration constant. Below we will show 
that this integration constant is to be interpreted as the 
mass of the central body. We will call the metric (|6"5]) 
together with condition (|T9"|) the McVittie metric. 

Clearly the system (|78|l is under-determining. This is 
expected since no equation of state has yet been imposed. 
The reason why we did not impose such a condition can 
now be easily inferred from (JT5J): whereas (|78b|) implies 
that g only depends on t, (|78c[) implies that p depends 
on t and r iff (a/a) =^ 0. Hence a non-trivial relation 
p = p(g) is simply incompatible with the assumptions 
made so far. The only possible ways to specify p are 
p = or q + p = 0. In the first case (|78c[) implies that 
a/a = if mo ^ (since then the second term on the 
right-hand side is r dependent, whereas the first is not, 
so that both must vanish separately), which corresponds 
to the exterior Schwarzschild solution, or a(t) oc i 2 / 3 if 



19 



mo = 0, which leads to the flat FLRW solution with 
dust. In the second case the fluid just acts like a cos- 
mological constant A = 8irg (using the equation of state 
g + p = in divT = it implies dp = and this, in 
turn, using again the equation of state, implies dg = 0) 
so that this case reduces to the Schwarzschild-de Sitter 
solution. To see this explicitly, notice first that (|78b|78c|l 
imply the constancy of H = a/ a = \/A/3 and hence 
one has a(t) = oq exp (t y/A/3) . With such a scale- 
factor the McVittie metric ([55]) with ([79")) turns into the 
Schwarzschild-de Sitter metric (|59|) in disguise. The ex- 
plicit formulae for the coordinate trans formation r e lating 
the two can be found in Section 5 of dRobertsonl. Il928h 
and also in Section 7 of ( Klioner and Soffel . l2005f h Fi- 
nally, note from (|78a[) that constancy of one of the func- 
tions m and a implies constancy of the other. In this case 
(|78bl78c[) imply g = p = 0, so that we are dealing with 
the exterior Schwarzschild spacctimc. 

A specific McVittie solution can be obtained by choos- 
ing a function a(t), corresponding to the scale function 
of the FLRW spacetime which the McVittie model is 
required to approach at spatial infinity, and the con- 
stant mo, corresponding to the 'central mass'. Rela- 
tions (|78b|78cp . and (|79|) are then used to determine g, 
p, and m, respectively. Clearly this 'poor man's way' 
to solve Einstein's equation holds the danger of arriving 
at unrealistic spacetime dependent relations between g 
and p. This must be kept in mind when proceeding in 
thi s fashion. For furthe r discussion of this point we refer 
to (|Nolanl . ll998lll999aD . 



1. Interpretation of the McVittie model 

In this section we discuss the interpretation of the 
McVittie model, its singularities, trapped regions, sym- 
metry properties, and also the motion of the matter. In 
doing this, we shall take care to isolate those properties 
which arc intrinsic to the ansatz (|68|) independent of the 
imposition of Einstein's equation. The analysis can then 
also be applied to all generalizations which maintain the 
ansatz (|68|) . Gene ralizations in this sense ha ve recently 
been discussed in ( Faraoni and J acquest I2007T I. on which 
we will comment at the end of this section. 

According to what has just been said, we wish to re- 
gard the McVittie solution as a candidate model for an 
isolated mass mo in an 'otherwise' flat FLRW universe 
with scale function a(t). As already emphasized in the 
introduction, this requires specific justification in view of 
the fact that simple superpositions of solutions are dis- 
allowed by the nonlinearities. A set of criteria for when 
a solution represents a localized mass immersed in a fiat 
FLRW b ackground ha ve been proposed and discussed in 
detail in (jNolanl . Il998l ) . The basic idea is to employ the 
MS energy (in a spherically symmetric context, where it 
is equivalent to the Hawking mass) in order to detect lo- 
calized sources of gravity. We will follow this approach 
and for this purpose we compute the Ricci and the Weyl 



part of the MS energy. 

This we now do for the class of metrics with- 
out at first making any use of Einstein's equation. The 
geometric definition of the MS energy in terms of the 
sectional curvature, together with formula (|Al0b[) for its 
Ricci part specialized to metrics with spatially isotropic 
Einstein tensor, implies that the Ricci part of the MS 
energy is given by 



E R = ±i? 3 Ein(e ,e ) 



(80) 



The Weyl part is then obtained as the difference between 
the full MS energy and (|80|) . We use the expression ([58]) 
for the former and write (dR, dR) = (e (R)) — (ei (R)) . 
The part involving eo(i?) equals ([80]) . due to the relation 
Ein(eo,eo) = 3(dR(eo) / R) 2 , which, e.g., follows from 
the comment below Eq. (|75|) an d (|D54jl (for vanishing 
shear) . The Weyl part of the MS energy is therefore given 
by (R/2)(l - (ei(i?)) 2 ). From {74]) we calculate e^R) 
and hence obtain for the Weyl part of the MS energy: 



Evj = am . 



(81) 



Now we invoke Einstein's equation with source ([70]) 
and four- velocity ([7T|) . Then the Ricci and Weyl contri- 
butions to the MS energy can be written in the following 
form, also taking into account (|79|) . 



(82a) 
(82b) 



47T o 

E R = — R 3 g, 
E w =m . 



Identifying the gravitational mass of the central object 
with the Weyl part of the MS energy, its constancy means 
that no energy is accreted from the ambient matter. As 
regards the Ricci part, note that the factor (47r/3)i? 3 
in (|82ap is smaller than the proper geometric volume 
within the sphere of areal radius R. This can be at- 
tributed to the gravitational binding energy that dimin- 
ishes the gravitational mass of a lump of matter below 
the value given by the proper space integral of T(eo, eo). 
This is shown in more detail in Appendix ID.4[ in par- 
ticular in the exact equation (|D47[) and its leading order 
approximation (|D49[) . 

It is also important to note that the central gravita- 
tional mass in McVittie's spacetime may be modeled by a 
shear-free per fect-fluid sta r of positive homogeneous en- 
ergy density (|Nolanl . Il993) . The matching is performed 
along a world-tube comoving with the cosmological fluid, 
across which the energy density jumps discontinuously. 
This means that the star's surface is comoving with the 
cosmological fluid and hence, in view of (|8"9"|) . that it ge- 
ometrically expands (or contracts). This feature, how- 
ever, should be merely seen as an artifact of the McVittie 
model (in which the relation ([55)) holds), rather than a 
general property of compact objects in any cosmological 
spacetimes. Positive pressure within the star see ms to be 
only possible if 2ad + a 2 < (see Eq. (3.27) in ( Nolanl . 
I1993D with a = exp(/3/2)), that is, for deceleration pa- 
rameters q > 1/2. 
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Next we comment on the singularity properties of the 
McVittie model. From (|78cp it is clear that, unless a/ a 
is constant (the Schwarzschild-de Sitter case) or m = 
(FLRW case), the pressure diverges at r = m/2 (that is 
at R = 2juq = Rs)- In fact, this corresponds to a gen- 
uine curvature singularity which is built into the McVit- 
tie ansatz (|68[) independently of any further assumption. 
To see that r = m/2 (corresponding to R = 2am = 2E W ) 
is a singularity it suffices to consider the scalar curvature 
(i.e. the Ricci scalar) of (|55|) . 



Seal = -12F 2 - 6-^, 



(83) 



whic h is readily comp uted from (|72[) . 
In ( Carrera and Giulinil . [2009) we show that this 



becomes singular in the limit r — > m/2, with the 
only exceptions being the following three special 
cases: (i) m = and a arbitrary (FLRW spacetimes), 
(ii) a and m constant (Schwarzschild spacetime), and 
(Hi) (am)' = and (a/a)' = (Schwarzschild-de Sitter 
spacetime). This means that, as long as we stick to the 
ansatz at r — m/2 there will always (with the only 
exceptions listed above) be a singularity in the Ricci 
part of the curvature and thus, assuming Einstein's 
equation is satisfied, also in the energy momentum 
tensor, irrespectively of the details of its underlying 
matter model. Hence any attempt to eliminate this 
singularity by maintaining the ansatz (|68[) and merely 
modifying the matter model is doomed to fail. In 
par ticular, this is true fo r the g eneralizations presented 
in ( Faraoni and Jacauesl [2007), contrary to what is 
claimed in that work. We also remark that it makes 
no sense to absorb the singular factors X/B in front of 
the time derivatives by wr iting (A/B)d/dt as en an d 
then argue, as was done in (Far aoni and Jacquesl l2007h . 
that this eliminates the singularity. The point is simply 
that then eo applied to any continuously differentiable 
function diverges as r — > m/2. Below we will show that 
this singularity lies within a trapped region. Turning 
back to the McVittie model, recall that in this case it is 
assumed that the fluid moves along the integral curves 
of d/dt, which become lightlike in the limit as r tends 
to m/2. Their acceleration is given by the gradient 
of the pressure, which necessarily diverges in the limit 
r — > m/2, as one explicitly sees from (|9T|). For a more 
det ailed study of the geometric singularity at r = m/2, 
see (|Nolaiill999a llbh. 

For spherically symmetric spacetimes the Weyl part of 
the curvature has only a single independent component 
(see (|D8a|) ) which, by its very definition, is ~2/R 3 times 
the Weyl part of the MS energy (see (|D36[) ). The square 
of the Weyl tensor for the ansatz (ffi5|) may then be con- 
veniently expressed as (see (|5T|) and (|D37I0 



(Weyl, Weyl) = 48 



(am) 2 
— R6~ 



(84) 



covered by our original coordinate system, for which r > 
m/2 (that is R > 2E W ). 

It is instructive to also determine the trapped regions 
of McVittie spacetime. We do this just using the McVit- 
tie ansatz (|68|) and making no further assumptions. Re- 
call that a spacelike two-sphere S is said to be trapped, 
marginally trapped, or untrapped if the product 9 + 8~ of 
the expansions (defined below Eq. (|D29[) ) for the ingo- 
ing and outgoing future-pointing null vector fields nor- 
mal to S is positive, zero, or negative. Taking S to be 
Sr, that is, a sphere of symmetry with areal radius R, 
it immediately follows from (|D31[) that Sr is trapped, 
marginally trapped, or untrapped iff (dR, dR) is posi- 
tive, zero, or negative, respectively. This corresponds to 
timelike, lightlike, or spacelike dR, or, equivalently, in 
in view of to 2E — R being positive, zero, or neg- 
ative, respectively. Using ([80)) together with (|72a[) , the 
MS energy for the McVittie ansatz can be written as 
E = £ w + R 3 F 2 /2, so that 



2E-R = F R -R + R s , 



(85) 



where Rs-=2E W denotes the 'generalized' Schwarzschild 
radius. Note that in general F depends itself on the ra- 
dial coordinate — except for the McVittie case, in which 
one has F = H =: 1/Rh (Rh denotes the Hub- 
ble radius). For computational simplicity we special- 
i ze in the following to the McVittie case, referring to 
(|Carrera and Giulinil . [20091) 'or the general case. Doing 
this, (155)) becomes a cubic polynomial in R which is pos- 
itive for R = and tends to ±oo for R — > ±oo. Hence it 
always has a negative zero (which does not interest us) 
and two positive zeros iff 



Rs/Rh < 2/3^3 re 0.38. 



(86) 



This shows that R = also corresponds to a genuine 
curvature singularity, though this is not part of the region 



This clearly corresponds to the physical relevant case 
where the Schwarzschild radius is much smaller than 
the Hubble radius. One zero lies in the vicinity of 
the Schwarzschild radius and one in the vicinity of the 
Hubble radius, corresponding to two marginally trapped 
spheres. The exact expressions for the zeros can be easily 
written down, but are not very illuminating. In leading 
order in the small parameter Rs/Rh, they are approxi- 
mated by 

Rx « R s (I + (Rs/Rh) 2 ) , (87a) 
R 2 a R H (1 - R s /2R h ) • (87b) 

From this one sees that for the McVittie model the ra- 
dius of the marginally trapped sphere of Schwarzschild 
spacetime (Rs) increases and that of the FLRW space- 
time (Rh) decreases. The first feature can be under- 
stood as an effect of the presence of cosmological matter, 
whereas the latter is an effect of the presence of a cen- 
tral mass abundance. All the spheres with R < Ri or 
R > i?2 are trapped and those with R\ < R < R 2 are 
untrapped. In particular, the singularity r = m/2, that 
is R = 2E W = Rs, lies within the inner trapped region. 
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Another aspect concerns the global behavior of the 
McVittie ansatz (1551) . We note that each hypersurface 
of constant time t is a complete Riemannian manifold, 
which, besides the rotational symmetry, admits a dis- 
crete isometry given in (r, 9, <p) coordinates by 



(r,6,<p) = ((m/2) 



2 r ~l 



(88) 



This corresponds to an inversion at the two-sphere r = 
m/2 and shows that the hypersurfaces of constant t can 
be thought of as two isometric asymptotically-flat pieces 
joined together at the totally geodesic (being a fixed- 
point set of an isometry and hence also minimal) two- 
sphere r = m/2. Except for the time-dependent factor 
m(t), this is analogous to the geometry of the t — const, 
slices in the Schwarzschild metric (the difference being 
that (|88j) does not extend to an isometry of the space- 
time metric unless m = 0). This means, in particu- 
lar, that the McVittie metric cannot literally be inter- 
preted as corresponding to a point particle sitting at 
r = (r = is at infinite metric distance) in a flat 
FLRW universe, just like the Schwarzschild metric does 
not correspond to a point particle sitting at r = in 
Minkowski space. Unfortunately, McVit tie seems to have 
interpreted his solution in this fashion ( McVittie) . Il933h 
which even until re cently gave rise to some confusion in 
the litera t ure (e .g. ( Ferraris et all fl996t iGautreavj , 1984 ; 
SussmarJ . Il988f )). A clarification was given by (jNolan . 
1999^ 

We now briefly discuss the basic properties of the mo- 
tion of cosmological matter for the McVittie model. We 
already mentioned that the vorticity and shear of the 
four-velocity u vanish identically. On the other hand, 
the expansion (divergence of it) is 



9 = 3H, 



(89) 



just as in the FLRW case (recall that here H :— a/a is 
defined as in the FLRW case, see (ff7|). In particular, 
the expansion of the cosmological fluid is homogeneous 
in space. The expression for the variation of the areal 
radius along the integral lines of u (that is the velocity 
of cosmological matter measured in terms of its proper 
time and the areal radius) is also just as in the FLRW 



u(R) = HR, 



(90) 



which is nothing but Hubble's law. The acceleration of 
u, which in contrast to the FLRW case does not vanish, 
is given by 



mo / 1 + m/2r 
^2 \J_ _ m /2r 



ei 



(91) 



In leading order in mo/R this corresponds to the 
acceleration of the observers moving along d / dT in 
Schwarzschild spacetime (see Eq. (| 132[) with A = 0). 

We conclude this subsection by commenting on the 
attempts to generalize the McVittie model. The first ob- 
vious generalization consists in allowing a non-vanishing 



cosmological constant. This is however trivial, since it is 
equivalent to the substitution g — ► g + qa and p — > p+pa 
in (|78|) . where qa '■= A/87T and pa '■= — A/8tt are, re- 
spectively, the energy-density and pressure associated to 
the c osmological constant A. In (jFaraoni and Jacauesl . 
l2007h . attempts to generalize the McVittie model have 
focused on keeping the ansatz (|68|) and relaxing the con- 
ditions on the matter in various ways. More precisely, 
they show that the McVittie case may be generalized to 
allow for radial fluid motions relative to the e observer 
field (that is relaxing condition (fTTj) ). provided one also 
allows for a non-vanishing radial heat flow. Both energy 
flows are necessary in order to get new solutions consis- 
tent with the ansatz (|68|) . even though the two radial 
energy flows do not cancel in the energy balance. As a 
result, the Weyl part of the MS energy will now change 
in time so that the new solutions correspond to inho- 
mogeneities of variable strength due to accretion or loss 
of energy. For further analysi s of these solutions we refer 
to ( Carrera and Giulinil.l2009fl . Another generalization of 
the McVitt ie model, this time away from the ansatz (|68[) . 
is given in (Sul tana and Dverl . l2005f) . It is constructed by 
applying a particular (time-dependent) conformal trans- 
formation ^to_th£_Schwa£zschilds_^cetime. As we will 
show in ( Carrera and Giulinil . [2009), the metric so ob- 
tained cannot be written in the McVittie form (1681) . 
contr ary to what is suggested in ( Faraoni and Jacauesl . 
2007). The corresponding energy-momentum tensor (ob- 
tained via Einstein's equation) can be interpreted as a 
sum of two contributions, one due to a perfect fluid and 
the other to a null fluid. 



2. Motion of a test particle in McVittie spacetime 

We are interested in the motion of a test particle (ide- 
alizing a planet or a spac ecraft) in McVittie's space- 
time. In ( McVittie] . [l933) it was concluded within a 
slow-motion and weak-field approximation that Keple- 
rian orbits do not expand as measured with the 'cos- 
molo gical geodesic radius' r* = a(t)r Late r ( Pachneii 
Il963l ) and (jNoerdlinger and Petrosiaril . Il97l[ ) argued for 
the presence of the acceleration term ([7]) proportional 
to d/a within this approximation scheme, hence arriving 
at (|12ap . In the following we shall show how to arrive at 
(I12ap from the exact geodesic equation of the McVittie 
metric by making clear the approxim ations involved. Re - 
lated recent discussions were given in (|Bolen et aLl . l200lh . 
where the effects of cosmological expansion on the peri- 
astron precession and eccentricity are discussed for con- 
stant Hubble parameter H := a/a. 

We will again work with the areal radius R. Note 
that for fixed t the map r i— > R(t, r) is 2-to-l and that 
R > 2mo, where R = 2mp corresponds to r = mo /2a. 
Hence we restrict the coordinate transformation ([74)) to 
the region r > mo /2a where it becomes a diffeomorphism 
onto the regio n R > 2mn. (T he region R < 2mo was in- 
vestigated in ( Nolanl . Il999b').) Reintroducing factors of 
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c, McVittie's metric assumes the (non-diagonal) form in 
the region R > 2m (i.e. r > m /2a(t)) 

g = (1 - 2/xCR) - h(t, R) 
2h(t,R) 



2 ^ c 2 dt 2 



V 1 - MR) 



c dt dR 



dR 2 



1 - 2fi(R) 



where we put 
ft(R) 



mo 
R 



hit, R) 



H(t)R 



(92) 



(93) 



with H(t) := (a/a)(t), as usual. 

The equations for a timelike geodesic (i.e. parametrized 
with respect to proper time), t i— > z M (r) with g(z, z) — 
c 2 , follows via variational principle from the Lagrangian 
C(z,z) = {l/2)g^(z)z^z v . Spherical symmetry im- 
plies conservation of angular momentum. Hence we may 
choose the particle's orbit to lie in the equatorial plane 
6 = 7r/2. The constant modulus of angular momentum 
is 



R 2 ip = L . 



(94) 



The remaining two equations are then coupled second- 
order ODEs for t(r) and R(t). However, we may re- 
place the first one by its first integral that results from 
g{z,z) = c 2 : 

(l - 2fj,(R) - h 2 {t,R)) c 2 i 2 

2h(t,R) ■ R 2 
zd R — 



y/\-2n{R)" l-2/x(i?) 
The remaining radial equation is given by 

R - (l-2f,(R)-h 2 (t,R))^ 



o (95) 
~{L/R) 2 =c 2 . V 1 



(1 - 2n(R)) i 2 



1/2 



mo c 

R^H(t)(l-2fx(R)) 1/ 

+ H(t) 2 (l-n(R)-h 2 (t,R)))i 2 (96) 



(n(R)-h 2 (t,R)) R' 
1 - 2fi(R) R 
2^(R)-h 2 {t,R)) 



v 7 ! - MR) 



cH{t) (R/c)t 



0. 



Recall that mo = GM/c 2 , where M is the mass of the 
central star in standard units. 

Equations ()95|96|) are exact. We are interested in or- 
bits of slow-motion (compared with the speed of light) in 
the region where 



R s <R<Rh 



(97) 



Recall that Rs and Rh are the Schwarzschild and the 
Hubble radius, respectively (see ([I]),©). The latter con- 
dition clearly covers all situations of practical applica- 
bility in the Solar System, since the Schwarzschild radius 



R s of the Sun is about 3 km = 2-10~ 8 AU and the Hubble 
radius R H is about 13.7 • 10 9 ly = 8.7 • 10 14 AU. 

The approximation now consists in considering small 
perturbations of Keplerian orbits. Let T be a typical time 
scale of the problem, like the period for closed orbits or 
else R/v with v a typical velocity. The expansion is then 
with respect to the following two parameters: 



c ~ V RJ 



(98a) 



corresponding to a slow-motion and weak-field approxi- 
mation, and 



HT. 



(98b) 



corresponding to the approximation for small ratios of 
characteristic-times to the age of the universe. In order to 
make the expression to be approximated dimensionless, 
we multiply (95]) by 1/c 2 and (96]) by T 2 /R. Then we 
expand the right hand sides in powers of the parameters 
(|9"5)) . using the fact that h := (HR/c) w e x e 2 . From this 
and (MJ) we obtain (T2"|) if we keep only terms to zero- 
order in £j and leading (i.e. quadratic) order in e 2 , where 
we also re-express R as function of t. Note that in this 
approximation the areal radius R is equal to the spatial 
geodesic distance on the t — const, hypersurfaces. 

We already mentioned that in the special case 
of constant H = d/a the McVittie solution turns 
into the Schwarzschild-de Sitter metric (|59[) . which 
also describes the spacetime inside the Einstein-Straus 
vacuoles in case of non-vanishing cosmological con- 
stant. Recently, exact expressions in terms of hyper- 
elliptic integrals for the integrated geodesic equa- 
tion in Schwarzschild-de Sitter spacetimes where de- 
rived in ( Hackmann and Lammerzah l. 2008a, b). More- 
over, a general discussion of Solar-System effects in 
Schwarzschild-de Sitter, like gravitational redshift, light 
deflection, time delay, perihelion precession, geodetic pre- 
ces sion, and effects on Doppl er tracking, has been given 
in ( Kagramanova et allJ20(M) . For example, it was found 
that a non- vanishing A could account for the anomalous 
Pioneer acceleration if its value was — 10~ 37 m -2 , which is 
minus 10 15 times the current most probable value. That 
value would also give rise to a perihelion precession four 
orders of magnitude larger than the accuracy to which 
this effect has been measured today. 

We conclude by commenting on the geodesic 
equation in the generalizations o f McV ittie's 
mo del given in (iFaraoni and Jacauesl [2007) and 
in (Sultana and Dveii [2005) . An essential feature which 
distinguish these solutions from the McVittie one, is 
that in the former the Weyl part of the MS energy 
E w = am is not a constant as for McVittie but varies in 
time, meaning that there is an accr etion of cosmologica l 
matt er by the inhomogeneity (see ( Carrera and Giulinil . 
2009)). In view of the fact that the combination 



m/r = A 2 E W /R « E w /R 



(99) 
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contained in the McVittie ansatz gives (minus) the 'New- 
tonian' part of the potential in the slow- motion and weak- 
field approximation (see Section IV. D.2jl . we deduce that 
in order to get the geodesic equation for the generalized 
McVittie models it suffices to substitute mo with E w in 
the equation of motion derived in Section IV.D.2I This 
means that the strength of the central attraction varies 
in time, leading to an in- or out-spiraling of the orbits if 
E w is increasing or decreasing, respectively. 



3. Exact condition for non-expanding circular orbits in McVittie 
spacetime 



Analogously to Section I1V.B1 where we ask whether 
there exist non-expanding circular orbits (i.e. of con- 
stant areal radius) of the electron-proton system in a 
FLRW spacetime, we now ask whether there exist non- 
expanding circular orbits of an (uncharged) test parti- 
cle around the central mass. The necessary and suffi- 
cient condition for this to happen follows from insert- 
ing R = const in the radial part of the geodesic equa- 
tion (|96[) and using the normalization condition (|95p of 
the four-velocity in order to eliminate t. In terms of the 
dimensionless quantities h(t) := RH(t)/c, I := L/Rc, 
and n := mo/R, the condition for the existence of non- 
expanding circular orbits can be given the following form: 

Rj_ ( M (l + 3Z a )-Z 2 -fe 2 ) 

(i + iVr=2M ' 1 j 

As for the electron-proton system in an FLRW space- 
time (see (|51l) ) this is a first-order autonomous ODE for 
h{t) and therefore the Hubble function. In the present 
case the ODE is even simpler since it has the elementary 
form h = p(h 2 ), where p is a polynomial of degree two 
with constant coefficients. From (jlOOl) . to leading order 
in the small quantities fx, I 2 , and h 2 , we get the same 
approximate ODE (|54|) and hence the same approximate 
solutions (|55p . Hence, the same conclusions as drawn for 
the electron-proton system in FLRW apply here. 

From (|100p it follows that stationary solutions hit) = 
const =: ho, corresponding to an exponentially-growing 
scale factor (|52|) (and hence leading to a Schwarzschild 
de Sitter spacetime), are those where ho satisfies 



l 2 + h 2 
(1 + 3Z 2 ) 



//.. 



(101) 



where we used that the first factor on the numerator of 
the right-hand side of (|100p is nonzero, as can be imme- 
diately inferred from the normalization condition (|95|) . 
Notice that for a vanishing Hubble parameter (that is for 
ho = 0) the above condition reduces to the third Kepler 
law in Schwarzschild spacetime, as expected. The effect 
of a non-vanishing Hubble parameter is again that we 
must provide the orbiting particle with a smaller angular 
velocity (smaller /) in order to keep it on a non-expanding 
circular orbit with the same radius. The largest radius 



at which in a McVittie spacetime with exponentially- 
expanding scale factor (that is a Schwarzschild-de Sitter 
spacetime) there is a non-expanding circular orbit follows 
from poip in the limit I — > 0. Then the condition reduces 
to h^ = n which, solving for R, gives (RsRh/^) 1 / 3 . This, 
exactly corresponds to the critical radius (|16p , taking into 
account that qo — — 1 for an exponentially-growing scale 
factor. 



VI. KINEMATICAL EFFECTS 

In this section we discuss the influence of cosmic ex- 
pansion upon measurements of relative distances, veloc- 
ities, and accelerations. These kinematical notions loose 
their a priori meaning in general spacetimes, in particu- 
lar in time-dependent ones. Hence it is of utmost impor- 
tance to carefully reconsider statements concerning such 
notions and their precise relations to locally observable 
quantities. 



A. Einstein- versus cosmological simultaneity 

Misidentifications in the notion of simultaneity can 
give rise to apparent anomalies in velocities and ac- 
cele ration. Such an effect ha s e.g. been suggested 
in (IRosales and Sanchez- GomezL Il998h and again in 
(|Rosalest 12002) to be able to account for the PA. Their 
argument says that in a spatially flat FLRW universe 
the mismatch between adapted cosmological coordinates 
on the one hand and radar coordinates on the other just 
amount to an apparent difference in radial acceleration of 
magnitude |T]). We agree on the existence and conceptual 
importance of such an effect but we disagree on the mag- 
nitude, which seems to have been grossly overestimated 
as we will show below. 

The cause of such effects lies in the way one actually 
measures spatial distances and determines the clock read- 
ings they are functions of (a trajectory is a 'distance' for 
each given 'time'). The point is this: equations of mo- 
tions give us, for example, simultaneous (with respect 
to cosmological time) spatial geodesic distances as func- 
tions of cosmological time. This is what we implicitly 
did in the Newtonian analysis. But, in fact, spacecraft 
ranging is done by exchanging electromagnetic signals. 
The notions of spatial distance and simultaneity thereby 
implicitly used are not the same as those we referred to 
above. Hence the analytical expression of the 'trajectory' 
so measured will be different. 

We first recall the local version of Einstein simultaneity 
in general spacetimes (./Vf,t/). We take ds — g^u dx^dx v 
to carry the unit of length so that dr = ds/c carries 
the unit of time. In general coordinates {x^} = {t, x 1 }, 
where x° = t denotes the timelike coordinate, the metric 
reads 



ds 2 = g^udx^ dx u = gttdt 2 + 2g t idt dx % 



gij dx l dx J . 

(102) 
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The observer at fixed spatial coordinates is given by the 
vector field (normalized to g(u,u) = c 2 ) 

u = c \\aiat\\- x a/at = -^=8/ at . (103) 

V9tt 



Consider the light cone with vertex p £ Ai; one has ds 2 = 
0, which allows to solve for dt in terms of the dx 1 (all 
functions g^ v are evaluated at p, unless noted otherwise) : 



dt 12 = ~^ dx l ± 
9tt 



9u9tj 9ij 



9tt 



^ ) dx i dx3 . (104) 

9tt, 



The plus sign corresponds to the future light-cone at p, 
the negative sign to the past light cone. An integral line 
of u in a neighborhood of p cuts the light cone in two 
points, q+ and If t p is the time assigned to p, then 
tq + = t p + dti and t q _ = t p + dt 2 . The coordinate-time 
separation between these two cuts is t q+ —t q _ = dt\ — dt2, 
corresponding to a proper time y/gtt{dti — dt 2 )/c for the 
observer u. This observer will associate a radar-distance 
dl* to the event p of c/2 times that proper time interval, 
that is: 



dli = h = 



9u9tj 
gtt 



9ij 



dx 1 dx 3 



(105) 



The event on the integral line of u that the observer 
will call Einstein-synchronous with p lies in the middle 
between q + and g_. Its time coordinate is in first-order 
approximation given by \{t q+ J rt q _) = t p + ^(dti+dt2) = 
t p + dt, where 

dt := Udtx + dh) = - 9 —dx l . (106) 
9tt 

This means the following: The integral lines of u 
are parametrized by the spatial coordinates {a; l }i=i,2,3- 
Given a point p, specified by the orbit-coordinates x p and 
the time-coordinate t p , we consider a neighboring orbit 
of u with orbit-coordinates x l p + dx 1 . The event on the 
latter which is Einstein synchronous with p has a time 
coordinate t p + dt, where dt is given by (|106[) . or equiva- 
lently 



9:=dt 



9u 
flu 



dx' = 0. 



(107) 



Using a differential geometric language we may say that 
Einstein simultaneity defines a distribution 9 = 0. 

The metric (I102[) can be written in terms of the radar- 
distance metric h (|105[) and the simultaneity 1-form as 
follows: 



ds 2 = g^ydx^ dx v = gtt 2 — h , 



(108) 



showing that the radar-distance is just the same as the 
Einstein-simultaneous distance. A curve 7 in M. inter- 
sects the flow lines of u perpendicularly iff 9(j) = 0, 
which is just the condition that neighboring clocks along 
7 are Einstein synchronized. 



We now apply the foregoing to isotropic cosmological 
metrics. In what follows we drop for simplicity the an- 
gular dimensions. Hence we consider metrics of the form 



ds 2 = c 2 dt 2 



a{t) 2 dr 2 



The comoving observer field, 

u = cd/dt , 



(109) 



(110) 



is geodesic and of expansion 3H. On a hypersurface of 
constant t the radial geodesic distance is given by a(t)r. 
Making this distance into a spatial coordinate, r», we 
consider the coordinate transformation 



t h- ► t* := t , v 1 > r% :— a(t)r . 
The field d/dt* is given by 

8/dU = 8/dt- Hrd/dr , 
to which the observer field, 

m* := cWd/auw- 1 d/dU 



(111) 



(112) 



(113) 



corresponds. In contrast to (|110[) . whose flow connects 
comoving points of constant coordinate r, the flow of 
(|112p connects points of constant geodesic distances, 
as measured in the surfaces of constant cosmological 
time. This could be called cosmologically instantaneous 
geodesic distance. It is now very important to realize that 
this notion of distance is not the same as the radar dis- 
tance that one determines by exchanging light signals in 
the usual (Einstcinian) way. Let us explain this in detail: 
From (|111|) we have adr = dr* — r*Hdt, where H := 
a I a (Hubble parameter). Rewriting the metric p09[) in 
terms of and yields 



ds 2 = c 2 (l - (Hr^/c) 2 ) dtl - drl + 2Hr* dtdr* 



c 2 {l-(Hr*/c) 2 }{ 



dt„ 



Hr,/c 2 



1 - {Hr*/c) 2 



dr* 



9t,t, 
drl 



1 - (Hr*/c) 2 ■ 



(114) 



Hence the differentials of radar-distance and time-lapse 
for Einstein-simultaneity are given by 



dL = 



dt* 



dr* 



y/\ - (Hr*/c) 2 ' 
Hr*/c 2 



1 - (Hr*/c) 2 



dr* 



(115a) 
(115b) 



Let the distinguished observer (us on earth) now move 
along the geodesic r» = 0. Integration of (|115p from = 
to some value r* then gives the radar distance l„ as well 
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as the time lapse At* as functions of the cosmologically 
simultaneous geodesic distance r*: 

/* = (c/H) arcsin(_ff r*/c) 

= r*{l + \{HrJcf + 0({HrJc) 3 )} (116a) 
At* = (1/2H) ln(l - (ff r*/c) 2 ) 

= (r*/c){-±(ffr*/c) + 0((ffr*/c) 2 )} (116b) 



observer: r» 



Combining both equations in (|116[) allows to express the 
time-lapse in terms of the radar-distance: 



At* = H- 1 ln(cos(in*/c)) 



(117) 



Now, suppose a satellite S moves on a worldline r* (t* ) 
in the neighborhood of our worldline r* = 0. Assume 
that we measure the distance to the satellite by radar 
coordinates. Then instead of the value r* we would use 
h and instead of the argument t* we would assign the 
time t* — At* which corresponds to the value of cosmo- 
logical time at that event on our worldline that is Einstein 
synchronous to the event (t*,r*); see Fig. [2] Hence we 
have 

l*(U) = ( C / H ) sin" 1 !^* + AQH/c} (118a) 



r* - ±(v/c){Hc){r*/cf 



(118b) 



where (I118b|l is (|118ap to leading order and all quantities 
are evaluated at t*. We set v = r*- 

To see what this entails we Taylor expand in t* around 
t* = (just a convenient choice): 



r*(t*) — r + v t* + \a Q t 2 + 
and insert in (fTT8b| . This leads to 



/»(t*) = r + u t* + 5O t* + 



where, 



f = r - (ffc) i(w /c)(r /c) z 
wo = «o - (#c) (v Q /c) 2 (r /c) 



(119) 



(120) 



(121a) 
(121b) 




satellite: r*(t») 



t* = const. 



FIG. 2 An observer moves on the geodesic worldline r, = 
and observes a satellite by exchanging electromagnetic signals 
(dashed line). Event po corresponds to the signal emission by 
the observer, p\ to its reflection by the satellite, and p2 to 
its re-absorption by the observer, t* denotes the cosmological 
time and the two curved lines correspond to hypersurfaces of 
constant t*. On the observer's worldline event A is defined to 
be simultaneous to p\ with respect to t* and A' is defined to lie 
half-way in proper time between po and p2 on the observer's 
worldline. In cosmological time A' is At* ahead of A. In cos- 
mological coordinates, the satellite's trajectory is represented 
by the a function t* i— + r,(t*), where r*(t«) denotes the proper 
geodesic distance in the hypersurface t* = const, between 
its intersection points with both worldlines. However, using 
radar coordinates, the observer takes A' to be simultaneous 
with pi and uses I, as measure for the satellite's simultane- 
ous distance. Since r* and i* are related by (|116a|l . it follows 
that the observer uses the function t* ?*(r*(t* — At*)) to 
characterize the satellite's trajectory, which leads to (|118|) . 



of magnitude! 17 Hence we conclude that, with respect to 
the PA, there is no significant kinematical effect result- 
ing from the distinct simultaneity structures inherent in 
radar and cosmological coordinates. 



5 = a - (He) {(v /c) 3 + (r /c)(v /c)(ao/c)} . (121c) 

These are, in quadratic approximation, the sought-after 
relations between the quantities measured via radar 
tracking (tilded) and the quantities which arise in the 
(improved) Newtonian equations of motion (not tilded). 

In particular, the last equation (|121cp shows that there 
is an apparent inward pointing acceleration, given by He 

times the (v/c) 3 + ■ ■ ■ term in curly brackets. As dis- 17 q ui . E a iTT7\ corre s ponds 
cussed in the introduction, He is indeed of the same or- 
der of magnitude as the PA, as was much emphasize d 
in (jRosales and Sanchez-Gomezl Il998t iRosalesl 120021) . 
However, in contrast to these authors, we also get the 
additional term in curly brackets, which in case of the 
Pioneer spacecraft suppresses the He term by 13 orders 



to Eq. 
From it 



(10) of 
the authors 



( Rosalcs and Sanchez-Gome 3. 1 1 99811. 

of jRosales and Sanchez-Gomej . I1998T) and llRosalesl . l2002t) 
immediately jump to the conclusion that there is "an effective 
residual acceleration directed toward the center of coordinates; 
its constant value is He". We were unable to follow this 
conclusion. Likew i se, we are unable to follow the conclusion in 
l lFahr and Siewertl. I200S?) . 
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B. Doppler tracking in cosmological spacetimes 

Doppler Tracking is a common method of tracking the 
position of vehicles in space. It involves measuring the 
Doppler shift of an electromagnetic signal sent from a 
spacecraft to a tracking station on Earth. This signal 
is either coming from an on-board oscillator or is coher- 
ently transponded by the vehicle in response to a signal 
received from the ground station. Here we focus on the 
second of these modes, which is more useful for navi- 
gation, partly because the returning signal is measured 
against the same frequency reference as that of the orig- 
inally transmitted signal and partly because the Earth- 
based frequency reference is also more stable than the 
oscillator on-board the spacecraft. 



1. Minkowski spacetime 

It is clear that this method will be fundamentally in- 
fluenced if performed within a time varying background 
geometry. Before elaborating on this, we consider the 
simple case of static Minkowski space. 



in 



V 



vehicle 



observer 




Po 



FIG. 3 Exchange of electromagnetic signals (represented by 
their rays at a slope of 45 degrees) between us and the space 
vehicle. Time runs vertically. 

In Fig. [3] we depicted two worldlines, one of the ob- 
server (straight vertically) and one for the space vehicle. 
A light signal is emitted by the observer at the event po, 
reflected by the vehicle at event pi, and finally received 
back by the observer at event p 2 - We choose a global 
Minkowski frame, that is global coordinates {x^} — 
{t,x 1 } with g(d/dx», d/dx v ) = diag(l, -1, -1, -1), in 
which the observer (for simplicity assumed to be iner- 
tial) is at rest at the origin of the spatial coordinates. If 



(3 := v/c denotes the radial velocity of the vehicle in units 
of c, the well known special-relativistic Doppler formula 
(applied twice) says that the ratio between the received 
and the emitted frequencies is 18 



LQ2(t 2 ) = 1 -/?(*! ) 



(122) 



Here to,ti, and t<i refer to the global Minkowski time at- 
tributed to po,Pi, and p%, respectively. In Doppler track- 
ing one is interested in the derivative of this ratio with 
respect to ti, which yields a measure for the velocity of 
the spacecraft. We assume u>q to be constant in time and 
note that, given the worldlines of the observer and the 
vehicle, t\ and to are uniquely determined by ti (since 
the events p\ and po are determined by p^). If r denotes 
the spatial radius coordinate we have t% — t\ = r(t\)/c. 
Differentiation with respect to ti leads to 



~d±2 



1 



1+J0(ti) 



and hence 



\ -3 



lW)= Mh){i+my 



(123) 



(124) 



This shows that — W2/2wo, namely (minus one-half) the 
derivative of the received to emitted frequency ratio with 
respect to the proper time of the receiving observer, gives 
the spacecraft's spatial acceleration up to corrections of 
order (3. Note that in view of note 1181 it would be inap- 
propriate to call these corrections 'special relativistic'. 

The final goal of this section is to derive the general- 
ization of (|124j) for a cosmological spacetime. For this we 
need two things: First, we need to know what is the gen- 
eralization of the concepts of spatial velocity and spatial 
acceleration in an arbitrary spacetime and, second, we 
need to know how electromagnetic signals propagate in 
an arbitrary spacetime. This is taken care of in the next 
paragraph. 



2. General setting 

In order to generalize the notions of spatial velocity 
and spatial acceleration to arbitrary spacetime one needs 



18 Note that in Special Relativity the Doppler Formula does, 
of course, not distinguish between moving emitter and mov- 
ing receiver. So l|122jl is obtained by squaring the frequency 
shift ^(1 - /3)/(l + /3), which is picked up once for the ra- 
tio uj fi(t r) / uji (ti) (receiver moving relative to the Minkowski 
frame) and once for W2(t2)/<*>i(ti) (emitter moving relative to 
Minkowski frame). Incidentally, exactly the same formula would 
result in non-relativistic physics if the observer is taken to be 
at rest with respect to the wave-guiding medium (e.g., the 
ether), which distinguishes the two states of relative motion. 
Indeed, in this case we have u\(t\)/u)o(to) = (1 — (3(t\)) and 
^2 (*2)A^i (ii ) = 1/(1 — /3(il)), whose product is again just H122I I, 
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to int roduce a fiducial reference ' observer-field': compare, 
e.g., (|Bini et aLl . ll995h and also (jCarreral [20091) ). An ob- 
server at the event p is a future pointing unit timelike 
vector in the tangent space T p {M.) of M. at p. An ob- 
server field is a field of future pointing unit timclikc vec- 
tors. Any observer u at p gives rise to an orthogonal split 
of the tangent space T p (A4) at p in a part parallel to u 
(the local time-axis) and a part orthogonal to it (the local 
rest space). Since u is not lightlike the two orthogonal 
subspaces are complementary, that is, together they span 
the whole tangent space and intersect only in the zero 
vector. The orthogonal projections of an arbitrary vec- 
tor X E T p (A4) onto these subspaces are, respectively, 
given by 

Q u {X):=g(X,u)u, (125a) 
P U (X) := X -g(X,u)u, (125b) 

which imply the decomposition identity X = Q U {X) + 
Pu{X). 

If two observers u and v are defined at the same point, 
the spatial velocity (over c) of v with respect to u is given 
by 



Pu(v) := 



P u (v) 

\Q u (v)\ 



v — g{v, uju 
g(v,u) 



(126) 



which is an element of the local rest space P U T(M). Its 
modulus is given by 



Pu(v) := \\Pu(v)\\ = ^l-l/g{u,vy. 



(127) 



Note that for the modulus we have f3 u (v) — (3 v (u), 
though the vectors (3 u (v) and (3 v (u) are linearly indepen- 
dent as they lie in P U T(A4) and P V T(A4), respectively. 
Note also that g(u, v) = 1/ \Jl — f3^(v) is just the ordi- 
nary 'gamma-factor'. Finally, if e e P U T(A4) is a unit 
vector, we define the spatial velocity of v in direction e 
w.r.t. it by 



/3» = -fl(e, /*») = - 



g(e,v) 
g(u,v) 



(128) 



The spatial acceleration of a worldline 7 w.r.t. a given 
observer field u is defined as the rate of change of the spa- 
tial velocity (3 u {j) within the local rest spaces P U T(A4) 
of u and with respect to the clocks moving along u. De- 
noting this acceleration (divided by c) with a, we have 



««(7) :=Vm.(7), 



(129) 



where we used the following covariant derivative for 
P U T(A4) -valued vector fields along 7: 



V? := \\QuM\~ 1 P u oV^oP u . 



(130) 



Here V-y denotes the ordinary (Levi-Civita) covariant 
derivative along 7. As an application one can, for ex- 
ample, rewrite the geodesic equation, V^,7 = 0, for a 



worldline 7 in terms of the spat i al qua ntities just intro- 
duced. One gets ( see (ICarreral 20091 ) or, in a slightly 
different notation, (jBini et all ll995f )T: 



'/3- u 



9u (/3k 



(131) 



where for better readability we omitted the arguments 
7 and 7 in the spatial acceleration and spatial velocity. 
Here a u := V u tt is the four-acceleration of the observer 
field it, U and ui u are, respectively, its shear-expansion 
and rotation tensors of rank (1, 1) (endomorphism), and 
Sp u :— P u + (3 U ® [3 U is a rank (1,1) tensor which, in 
a slow-motion approximation (neglecting quadratic and 
higher terms in j3), reduces to the identity on the lo- 
cal rest space of u. Equation (|13ip should be seen as 
a local version of Newton's equation. For example, in 
Schwarzschild-de Sitter spacetime (|59|) . taking u to be 
proportional to the timelike Killing field d/dT, one has 
9 U = = (because of the Killing equation and spher- 
ical symmetry, respectively) and 



a u = V u u — 



1 



m A 



(132) 



where eji denotes the the normalized radial vector field 
d/dR (we use here the coordinates and the notation 
of (I59|)). Hence, in slow-motion and weak-field approxi- 
mation (that is keeping only linear terms in /3, m/ R, and 
Ai? 2 ), the geodesic equation of motion in the form (|13ip 
reduces to 



m 
R~ 2 



, 2 11 : e,. . 



(133) 



which just gives the 'improved' Newtonian equation for 
geodesic motions in Schwarzschild-de Sitter spacetime. 
It has the same form as the improved Newtonian equa- 
tion studied in Section HTT1 

We turn now to electromagnetic signals and restrict 
our attention to monochromatic waves in the geometric- 
optics approximation (i.e. for wave-lengths negligibly 
small w.r.t. a typical radius of curvature of the space- 
time and w.r.t. a typical length over which amplitude, 
polarization, and frequency vary). In this approxima- 
tion an electromagnetic signal propagates on a light- 
like geodesic along which the wave-vector, fe, is tangent, 
future-pointing, and parallelly transported. Recall that 
fc is so normalized that the frequency measured by an 
observer, say it, is 



u u {k) := g(u, k) 



(134) 



Given a wave- vector fe and two observers u, v at the 
same spacetime point, their observed frequencies are thus 
uj v (k) — g(v, fc) and uj u {k) = g(u, fc), and their ratio is 
given by 



oj v (k) g(Q u v + P u v,k) 



uj u (k) 



g{u,k) 



g(u,v)[l-0t(v)] 



(135) 
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Here the spacelike unit vector fc := ||P„(fc)| 1 P u (k) de- 
fines the direction of k in the local rest space of u. In 
deriving ([135]) we used (TT28| and ||P«(fc)|| = g(u 7 k) 
to write g(v,P u (k)) = —g(v,u)g(u,k){3*(v). Equa- 
tion (|135|) is the general form of the Doppler formula. 

Let now u be an observer field along one integral line of 
which the distinguished observer is moving. The world- 
line of the vehicle is denoted by 7. The domain of the 
field u is assumed to include a neighborhood of 7. The 
wave- vector fc emitted at p suffers three changes: 

1. propagation from po to p\. feo — ► fci; 

2. reflection at p\. k\ — > k[; 

3. propagation from pi to p 2 : k[ — > fc 2 . 

We are interested in the ratio of the received to the emit- 
ted frequency: 



f2 _ g(u 2 ,k 2 ) 

uj g(u Q ,k Q ) 



w 2 



(136) 



What happens at reflection (the second process: fci — > 
k[)? Well, with respect to the spacecraft moving along 7 
with four-velocity v = 7, the wave vector fci at p\ splits 
according to 



fci = Qi(ki) + Py(ki) 



(137) 



A corner-cube reflector transported along 7 will reverse 
P-j(ki) while keeping Q-y(fei) intact (here we neglect a 
possible transponder shift which is irrelevant for our dis- 
cussion): 

fci 1 — ► k[ = Qy(ki) - Py(fci) = 2Q^(fci) - fci . (138) 

Hence u>\ := u> u (ki) = g(ui,k\) and lu[ := u> u (k[) = 
g(ui,k'i), the in- and out-going frequencies measured by 
the observer u at p\ , are related by 



-i = 2 



g{u,j)g{^f,k) 
g(u,k)\ Pl 



1-/52(7)U 



(139) 

where in the last step we just used (|135[) to rewrite the 
ratio g(j,k)/g(u,k). This accounts for the middle ratio 
on the right-hand side of (|136p . 

To account for the other two ratios in (|136[) . one uses 
the laws of geometric optics in (curved) spacetimc to re- 
late ll>q = g(uo, feo) (at p ) and lu 2 = g{u 2 , k 2 ) (at p 2 ) to 
kinematical quantities of 7 at pi. 19 For example, if u is 
a Killing field (like u = d/dt in Special Relativity), we 
haveg(w ,fco) = fl(«i,fci) and g(u 2 ,k 2 ) =g{u 1 ,k' 1 ), so 
that 



^2 _ 2 1-/3^7)U 1 
^0 1-^(7)U 



(140) 



As a trivial application, this includes the generalized form 
of (|122p . the latter corresponding to purely radial motion. 



19 In general spacetimes without timelike conformal Killing fields 
these quotients will also explicitly depend on time. 



3. FLRW spacetimes 

In standard cosmological spacetimes (FLRW), u = 
d/dt is not Killing, though X = a(t)d/dt is confor- 
mally Killing (Lxg = 2dg). One now has aog(uo, fco) = 
ai<7(«i,fci) and a 2 g{u 2 ,k 2 ) — aig(iii, fc'i), so that in- 
stead of (|140[) one gets 



« 2 1" i-^(t)I Pi 



0J 



- 1 



(141) 



We now want to relate the ^-derivative of (| 141 [) to 
the acceleration of 7. In order to calculate the derivative 
co 2 (t 2 )/uJo(to) we need to know the derivatives dt\/dt 2 
and dto/dt 2 . Restricting to the flat FLRW case for sim- 
plicity, they follow from the law of null propagation: 

f t2 dt 

/ — T-r = — — / dr, (142a) 



dt 

/o(f 2 ) a (*) 




ri(ti(t 2 )) 



dr— / dr 
'n(ti(t 2 )) J 



(142b) 

Differentiation with respect to t 2 yields, respectively 
dti a{ti) 



dt 2 

dtp 
dt 2 



(l + /t(7)U 
1-^(7)U 



a(t 2 ) 
a{to) 
a(t 2 ) 1 + ^( 7 )U 



(143a) 



(143b) 



The exact formula for the i2~derivative of the 
frequency-shift rate can now be computed. One obtains 




(144) 



where we suppressed the argument 7 and index u at 
(3 for better readability. This formula provides an ex- 
act relation between the time derivative of the observ- 
able frequency shift (defined 'here') and the kinemati- 
cal quantities of the vehicle (defined 'there'), provided 
the scale function a(t) is known. For purely radial mo- 
tion V^fc = and we obtain the simpler expression (now 



writing a for a k ) 
1 u 2 (t 2 ) 



a ai 



2 w (to) 



02 
a 1 



ao 
ai 



1-/3 
1 + /3 



l + (3 



(145) 
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In order to consistently approximate this expression in 
terms of small quantities j3 and HAt, where At := 
(*2 — to)/2, we think of (|145j) as being multiplied with 
At and regard aAt as being of order [3. Then, keeping 
only quadratic terms in [3, linear terms in H At where 
At := (£2 — io)/2, and also mixed terms (3 HAt, we get 

_ 1 UJ2^2l ^ _ _ 3FA A + _ (M6) 

2 Wo (to J 

Hence we see that in this approximation there are two 
modifications, besides the — 3/3-term already familiar 
from (|124|) , due to cosmic expansion: First, there is an 
additional contribution —3HAt acting in the same way 
as the — 3/3-term. It can also be interpreted in the same 
fashion, as its corresponds to the velocity (over c) of H At 
that a comoving systems picks up during the time the 
signal went from the observer to the vehicle. Second, 
there is a constant contribution H/3 to accclcration/c, 
i.e. Hc(3 to acceleration, in a direction parallel to the ra- 
dial velocity (i.e. outward pointing if the vehicle recedes 
from the observer). Hence it acts opposite to the PA 
and is smaller in modulus by a factor of (3. Applied to 
the Pioneer spacecrafts, the H At-term amounts to a tiny 
'anomalous' acceleration of Aa/a < 10~ 12 , the Hc(3-teim 
to Aa/a < 10" 7 . 

A final point must be made regarding the choice of the 
reference observer-field on which the kinematic quantities 
related to the spacecraft (spatial velocity and spatial ac- 
celeration) and the electromagnetic signal (frequency and 
spatial propagation direction) crucially depend. In the 
Minkowskian case the reference field was just u = d/dt, 
which is inertial, that is, geodesic and of vanishing ro- 
tation, shear, and expansion. It is clear that in a gen- 
eral spacetime such observer fields do not exist and there 
is no natural choice to replace them. However, in the 
case of spherical symmetry there is, in fact, a distin- 
guished observer field, namely that one whose orbits lie 
within the timelike hypersurfaces of constant areal ra- 
dius and there run perpendicular to the orbits of the 
rotation group. This clearly defines a non-rotating and 
'non-expanding' (w.r.t. the areal radius) reference field. 
It is the normalization of the so-called Kodama vector 
field, which we discuss in detail in Appendix ID. 31 In a 
FLRW spacetime it is just given by (|113p . Notice that 
in the present case, where the hypersurfaces of constant 
cosmological time t are flat, the areal radius corresponds 
also to the proper distance. Hence the integral curves 
of it* intersect the hypersurface of constant cosmological 
time at constant spatial geodesic distance. More pre- 
cisely, the expansion and the shear scalar of it* are given 
by d* = R 2 HH/{1 - {RHff 2 and cr* = -0*/3, re- 
spectively, showing that they are of order H 3 which we 
neglect. In passing we remark that the expansion and 
shear of u» exactly vanish for the de Sitter case, whose 
metric in 'static' coordinates is given by (|59|) for m = 0. 
In this case d/dt* = d/dT, that is, it, is proportional to 
the timelike Killing vector field d/dT; see (J5HJ). Coming 
back to the general FLRW case, the acceleration of it* 



is given by a u , = (~Rd/a + R 3 H 4 )/(l - {RH) 2 f/ 2 e*, 
where e* is the unit vector field orthogonal to it* and 
to the two-sphere, pointing in positive radial direction. 
Hence, in the slow-motion and weak-field approximation 
of Eq. (|146[) . but keeping also quadratic terms in H , the 
geodesic equation in the form (|13ip w.r.t. the observer 
field it* reads as 

au,(7) ~ (^ r *) e * 01 ■ ( 14? ) 

This is just an alternative derivation of the acceleration 
term ([7]). We point out that had we we taken (jllpp as 
observer field we would have arrived at the equation of 
motion a u (j) w —Hf3 u (j) instead of (|147|) . that is, no 
acceleration term (JTJ) would have resulted. 

In the approximation within which (| 146|) is derived 
this equation remains valid if the quantities in it are re- 
interpreted so as to refer to it* instead of u. Hence 
we may sum up the situation by saying that equa- 
tions (|146[) and (|147[) give, respectively, the two-way 
Doppler-tracking formula and the 'Newtonian' equation 
in a FLRW spacetime within the mentioned approxima- 
tion. 



4. McVittie spacetime 

The same analysis can be generalized from the spatially 
flat FLRW spacetime ((30} to the spatially flat McVittie 
spacetime (fTiS)) . Here the observer moves along d/dt, 
which is not geodesic. The coordinate t does now not 
measure proper time, denoted by r, along the observer's 
worldline. The result corresponding to (|146p can now be 
stated as follows: 

- - ^4 ~ a(l - 3/3 - 3At(H - m c/R 2 )) + Hp . 
2 w (to) 

(148) 

Here R denotes the areal radius of the observer during the 
measurement. Note that even though it changes along 
the observer's worldline according to (f9TJ)) . we do not need 
to account for the corresponding change in Atitiqc/ R 2 
of (— 2ATrriQc/ R 2 )(H At) which is of subleading order. 
The additional term in (|148|) has a straightforward inter- 
pretation in terms of the acceleration that the observer 
necessarily experiences while keeping a constant radius 
R away from the central inhomogeneity. 

As for the FLRW case, we chose the observer field to 
which we refer the spatial quantities to be proportional 
to the Kodama vector field (along which the areal ra- 
dius is constant). Putting r*(t, r) :— A 2 (t,r)a(t)r and 
t*(t, r) :— t, a short computation shows that the vector 
field d/dt* is again given by (|112[) . In the slow-motion 
and weak-field approximation used in Section [V.D.21 the 
geodesic equation in the form (|13ip w.r.t. the observer 
field it* reads 

au, (7) ~ (-'** — ) e* o 7 , (149) 
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where again e* denotes the unit outward-pointing vec- 
tor field orthogonal to it, and to the two-spheres of 
symmetry. This is an alternative derivation of the im- 
proved Newtonian equation for the McVittie spacetime 
carried out in Section IV. D. 21 Notice that, again, within 
the approximations used, relation (|148|) remains valid if 
one refers the quantities to it* instead of to u. There- 
fore (|148p and (|149p give the two-way Doppler-tracking 
formula and the improved 'Newtonian' equation for the 
McVittie spacetime within the mentioned approximation. 

In the special case of purely radial motion, insertion 
of (|149[) into (|148|) leads to a formula predicting the two- 
way Doppler-shift rate in linear order in HAt and Too /V* , 
and quadratic order in (3 Ut (j): 



1^2(72) 
2 w (to) 



toq 



(l - 3/3*) 



(150) 



Hence there are two corrections to the Newtonian con- 
tribution. One is proportional to H and stems from the 
cosmological expansion, the other, already familiar from 
the special- relativistic treatment (I124[) , is independent of 
H and merely due to the finiteness of the propagation 
speed of light (recall note [T8|) . Their ratio is (up to a 
factor given by the square of the ratio of r* to the 
geometric mean of the Schwarzschild radius too and the 
Hubble radius c/H. The latter is of the order of 10 23 km, 
so that its geometric mean with a Schwarzschild radius 
of one kilometer is approximately given by 2400 astro- 
nomical units. The ratio of the effects is therefore of the 
order 10 -7 . Hence the cosmological contribution is negli- 
gible for any application in the Solar System as compared 
to the 3/3-correction. For the Pioneer 10 & 11 space- 
crafts we have a radial velocity of about 12 Km/s. This 
amounts to a 3/3-correction of magnitude 4 • 10 -5 times 
the Newtonian gravitational acceleration, in an outward- 
pointing direction. This is indeed of the same order of 
magnitude as the PA but directed oppositely. 



VII. SUMMARY AND OUTLOOK 

We think it is fair to say that there are no theoret- 
ical hints that point towards a dynamical influence of 
cosmological expansion comparable in size to, say, that 
of the anomalous acceleration of the Pioneer spacecrafts. 
There seems to be no controversy over this point, though 
for completeness it should b e mentioned that there ex- 
ist speculations (|Palld . |2005[ ) according to which it might 
become relevant for future missions. But such specu- 
lations are often based on models which are not eas- 
ily related to the intended ph ysical situation, like that 
of Gautreau (jGautreaul . 119841 ). Rather, as the (a/a)- 
improved Newtonian analysis in Section lHTl together with 
its justification given in the subsequent Sections shows, 
there is no genuine relativistic effect coming from cos- 
mological expansion at the levels of precision envisaged 
here. 



On the other hand, as regards kinematical effects, the 
situation is less unanimous. It is very important to un- 
ambiguously understand what is meant by 'mapping out 
a trajectory', i.e. how to assign 'times' and 'distances'. 
Eventually we compare a functional relation between 
'distance' and 'time' with observed data. That relation is 
obtained by solving some equations of motion and it has 
to be carefully checked whether the methods by which the 
tracking data are obtained match the interpretation of 
the coordinates in which the analytical problem is solved. 
In our way of speaking, dynamical effects really influence 
the worldline of the object in question whereas kinemat- 
ical effects change the way in which one and the same 
worldline is mapped out from another worldline repre- 
senting the observer. Here we have derived exact results 
concerning the influence of cosmic expansion on this map- 
ping procedure, which allow to reliably estimate upper 
bounds on their magnitude. They turn out to be too 
small to be of any relevance in current satellite track- 
ings, which is in accord with naive expectation but in 
contrast to some statements found in the literature. 

At this point it is useful to recall once more the 
general philosophy behind such statements: From the 
Einstein-Straus solution it is clear that local overden- 
sities inhibit cosmic expansion, or at least that part of 
it which is not due to a cosmological constant. Also, 
as already mentioned before, the effect of anisotropics is 
also to diminish the effect of global expansion (see, e.g., 
( Dominguez and Gaitd . l200lh ) . Hence calculating such 
an effect in simple models like the improved Newtonian 
equation discussed in Section HI.AI (backed up by the var- 
ious justifications we discussed in detail) clearly means to 
overestimate the impact of cosmic expansion in a realistic 
situation, where the single overdensity (e.g. representing 
the Sun) is surrounded by more overdense structures (the 
Solar-System environment, the Galaxy, etc.) with less 
symmetry. If this overestimation gives an already in- 
significant upper bound for the envisaged effect, we can 
conclude that it becomes even more insignificant in more 
realistic models. 

Satellite navigation is clearly not the only potential 
source of interest in the question of how local inho- 
mogeneities affect cosmological expansion. Many pre- 
dictions concerning cosmological data rely on computa- 
tions within the framework of the standard homogeneous 
and isotropic models, without properly estimating the 
possible effects of local inhomogeneities. Such an esti- 
mation would ideally be based on an exact inhomoge- 
neous solution to Einstein's equations, or at least a fully 
controlled approximation to such a solution. The dy- 
namical and kinematical impact of local inhomogeneities 
might essentially influence our interpretation of cosmo- 
logical observations. As an example we mention re- 
cent serious efforts to interpret the same data that are 
usually taken to prove the existence of a positive cos- 
mological co nstant A in a c ontext w i th rea l istic inho- 
moge neities dBuchertl l200d: iRasanenl 120061: IWiltshird . 
l2007t) . i.e. taking into account that cosmological param- 
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eters are dresse d ( Buchert and Carforal . 12003^. See also 
(iBuchertl . l2008h br a recent review. One might spec- 
ulate that the measured A can eventually be fully re- 
du ced to the action of i nhomogeneities, as suggested in 
in ( Wiltshire) . 12007 , 20081). Fo r an earlier advance in this 
direction, see (|Celerierll200f)h . 
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APPENDIX A: Notation, conventions, and generalities 

A model for spacetime consists of a tuple (A4,g), where 
M is a four-dimensional manifold and g a Lorentzian 
metric whose signature we take to be (+, — , — , — ), i.e. we 
use the 'mostly minus' convention. Throughout we de- 
note geometric objects, like tensor fields and covari- 
ant derivative operators, by bold-faced letters or words. 
The unique metric preserving and torsion-free covariant 
derivative associated with g will be denoted by V and 
the covariant derivative in the direction of a vector X 
by Vx- For a smooth tensor field T on M its covariant 
derivative VT defines a linear map, X 1— > V^T, from 
the tangent space to the tensor space at each point of 
Ai where T is defined. Since VT is again a tensor field 
(of rank (p,q + 1) if the rank of T was (p,q)) we can 
form VVT := V(VT). Note that (VVT)(X,Y) = 
VxVyT — Vv x yT. For a scalar function / on M. we 
have V/ = df, the ordinary exterior differential, and 
W/ = Hess(/), the Hessian of /. The metric g al- 
lows to uniquely associate to any vector X a linear form 
X := g(X, ■ ), called the dual (with respect to g) of X. 
The inverse of this map will be denoted by an overline. 
The gradient of a function / is then grad/ = df. The 
metricity of V implies that the latter commutes with 
the maps j_ and ~ e.g. it holds: VxH = V xY ■ The 
scalar product induced by g on the tensor bundle will be 
denoted by (-,•). 

Associated with any two linearly independent vectors 
X,Y at a point p € M. is a curvature endomorphism, 
R{X, Y), of the tangent space at p: 



R(X,Y)Z = (VVZ)(X, Y) - (VVZ)(Y,X) 



(Al) 



= VxVyZ - VyVxZ - V[ X ,Y]Z. 

The Riemann- or curvature tensor, Riem, is then defined 



Riem(W, Z, X, Y) := g(W, R(X, Y)Z) 



(A2) 



It is antisymmetric under the exchange X <-> Y or W 
Z and symmetric under the slotwise exchange of pairs 
(W, Z) <-*■ (X,Y). Moreover, the antisymmctrization 
over any three slots vanishes (first Bianchi identity). The 
Ricci tensor, Ric, is defined by the trace of the following 
endomorphism 



Ric(r, Z) := tr(X ^ Riem(X, Y)Z) , 



(A3) 



which is symmetric under exchange Y <-* Z . The scalar 
curvature is defined by taking the trace of Ric, also called 
the Ricci scalar, with respect to g (since Ric is not an 
endomorphism, we need the metric to define its trace) 



Seal = tr g (Ric) . 



(A4) 



Finally, the Einstein tensor is the following combination 
of Ric and Seal: 



Ein := Ric — i Seal g . 



(A5) 



Associated to any spacelike or timelike two- 
dimensional plane 77 in the tangent space at p e Ai is 
the sectional curvature. Its geometric interpretation is 
just that of the ordinary Gaussian curvature at p of the 
two-dimensional surface in A4 that is spanned by the 
geodesic curves through p tangent to 77. In terms of 
Riem it reads 



Riem(X, Y,X,Y) 
Q(X,Y) 



(A6) 



where X, Y are any two linear independent vectors in 77 
and 



Q(X, Y) : = g(X, X)g(Y, Y) - g(X, Y) 
= (g&g)(X,Y,X,Y). 



(A7) 



Note that \Q(X, Y)\ gives the square of the area of the 
parallelogram spanned by X and Y which is nonzero 
iff the considered plane is spacelike or timelike (non- 
degenerate). 

In (|A7|) we introduced the product 0, which is called 
the Kulkarni-Nomizu product. It is a symmetric bilinear 
map from the space of symmetric (0, 2) tensors to the 
space of (0, 4) tensors with the same algebraic symmetries 
as Riem. Its general definition is as follows: 

(aQb)(W,Z,X,Y) := 

l -(a{W,X)b{Z,Y)-a{W,Y)b{Z,X) (A8) 

+b(W, X)a(Z, Y) - b(W, Y)a(Z, X) 

This can be used to conveniently write down the g- 
orthogonal decomposition of the curvature tensor into 
the Ricci- and the Weyl part: 



Riem = Ricci + Weyl 



(A9) 
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In four spacetime dimensions one has 

Ricci : = (Ric-| Seal g) g (AlOa) 
= (Ein-±tr g (Ein)g) g . (AlOb) 

Inserting this into (|A9|) gives the definition of Weyl. The 
definition is such that the Ricci part is g-orthogonal to 
the Weyl part and that the latter is totally trace free. 
Hence the Ricci and the Weyl part each contribute 10 
independent components to the 20 independent compo- 
nents of Riem. The Ricci part may be further decom- 
posed according to the decomposition of Ric into its 
trace and a trace-free part, but this refinement will not 
be needed here. 

Einstein's equation now express the local determina- 
tion of the Ricci part of the curvature in terms of the 
energy-momentum distribution of matter, the latter be- 
ing encoded in the energy-momentum tensor T of the 
matter. In units where Newton's constant G and the 
velocity of light c equal one 20 , Einstein's equation reads 

Ein = 87rT. (All) 

Here we did not write down explicitly a cosmological 
term, which can always be thought of as extra contri- 
bution to T of the form gA/8n. Now, assuming that g 
satisfies Einstein's equation, the Ricci part of the Rie- 
mann tensor is given in terms of T by 

Ricci = 8tt(T- \tr g (T)g) Qg. (A12) 

APPENDIX B: Proof of Theorem [1] 

In this section we prove Theorem [TJ namely the equiv- 
alence, in the spherically-symmetric case, of the SSJC 
with the Darmois junction conditions. 

Proof. The proof essentially consists in writing down the 
induced metric and extrinsic curvature for a (non-null) 
spherically symmetric hypersurface in a spherically sym- 
metric spacetime. This is most easily done by introducing 
an adapted orthonormal frame. 

We first consider the case where r is timelike, hence 
7 = 7r(_r) is a timelike curve in B. The following con- 
struction shall be carried out in both spacetimes. One 
defines v as in the SSJC, hence as the (unique up to a 
sign) spherically symmetric, unit vector field on r or- 
thogonal to n. That is v, seen as a vector field on B, 
is tangent to 7. Since n is spacelike, v is timelike. The 
ambient metric can be then written as 

g = v_ v — n n — R 2 g s 2 , (Bl) 



so that the induced metric (compare Appendix [C] and 
(Ol) on r is 

g r = v® v — R 2 gs 2 ■ (B2) 

In view of (jC3[) note that here s(n) — — 1. For the ex- 
trinsic curvature (jC6[) , using (|D3|) and the fact that v is 
spherically symmetric and hence tangent to B, one has 
the decomposition 

Kp = -g(n, V v v) v v - RdR(n)g S 2 . (B3) 

Now, from expressions (|B2|) and (|B3[) it follows that 
the DJC, and hence the continuity of gp and Kp, are 
equivalent to the continuity of the following four func- 
tions: 

(a) the arc-length of 7, 

(b) R, 

(c) dR(n), and 

(d) g(n,V v v). 

The statement of the theorem will now follow from the 
following expression of the MS energy (|58[) : 

E=^(l + (dR(v)) 2 -(dR(n)) 2 ) . (B4) 

Simply note that if R is continuous through _T (recall 
the definition of this concept below the definition of DJC 
in Section IV. Ap the same holds for its derivatives tan- 
gent to r. In particular, dR(v) is continuous through r 
and hence we may substitute dR(n) by the MS energy 
in the above list (a)-(d). This completes the proof for 
timelike r. 

In the case of spacelike r the unit normal n is timclike 
and v is chosen as the unique (up to a sign) spherically 
symmetric unit vector field on r orthonormal to n. Then 
v is a spacelike 'radial' unit vector field orthogonal to 
the S'0(3)-orbits. The proof now proceeds analogously 
to the timelike case. We merely list the expressions for 
the ambient metric 

g = n®n-»®» - R 2 g s 2 1 

the induced metric 

g r = v v + R 2 g S 2 , (B5) 

the extrinsic curvature 

K r =g(n, V v v)v®v + RdR(n)g S 2 , (B6) 

and the MS energy 

E=*{l + (dR(n)f-(dR(v)f) , (B7) 



Otherwise the factor 87r on the right-hand side of 1A11^ should 
be replaced with SttG/c 4 . 



and conclude exactly as in the timelike case. 



□ 
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APPENDIX C: Submanifolds 

In a Lorentzian manifold ( M , g) endowed with Levi- 
Civita connection V consider a smooth submanifold r 
of co-dimension one and normal vector field n. We as- 
sume r to be non-null, that is, either spacelike (then n 
is timelike) or timelike (then n spacelike). Then r in- 
herits from the ambient manifold M. a (non-degenerate) 
metric and a connection in a natural way. We introduce 
the orthogonal projectors 



Qn := e(n) n ® n 
P n := id-Q n , 



(Cla) 
(Clb) 



where e(n) denotes the indicator, defined for any non- 
null vector by 



g(X,X) _ J +1 if X is timelike, 
\g(X,X)\ I — 1 if X is spacelike. 



(C2) 



The induced metric on f (also called first fundamental 
form) is given by 



g r ■= -e{n)P n g , 



(C3) 



where the sign is just in order to get a positive definite 
metric in the case where r is spacelike. Given two vector 
fields X,Y tangent to F, so that Q n X = Q n Y = 0, 
one may decompose the covariant derivative of Y with 
respect to X into its orthogonal components 

V X Y = P n (V x Y) + Q n (V x Y) 

= r V x Y + K r (X,Y)n, (C4) 



where 

r V x Y :=P n {V x Y) 

is the induced connection on J 1 and 

K r (X,Y) :=e(n)g(V x Y,n) 
= -e(n) g(V x n,Y) 



(C5) 



(C6) 



is the extrinsic curvature of _T in (also called the sec- 
ond fundamental form). The second equality sign in (|C6j) 
is an immediate consequence of the metricity of V and 
the fact that X and Y are orthogonal to n. With this 
alternative expression for the extrinsic curvature one has 
K i = —e(n) P n Vn and hence 21 



K r = -e(n)S(P n Vn) 



(C7) 



Since n is hypersurface orthogonal (by definition) we 
have yl(P n Vn) = 0. Hence, the extrinsic curvature is a 



21 Here and below S and A denote the projection operators of full 
symmetrization and full antisymmetrization, respectively. 



symmetric (0, 2)-tensor field. We recall also that the in- 
duced connection is the Levi-Civita connection of (r, gr), 
as one may easily check. 

The full relations between the curvature of M. and 
those (i ntrinsic and extrinsic) of T can be found, 
e.g., in (jGiulinil . Il998l ). Here we are only interested in 
the 'Einstein part' of the curvature. One gets: 

Ein(n,n) = - (-e(n) r Scal +(trK) 2 - ||K|| 2 ) (C8a) 
Ein(n, P n ) = -e{n) div r (K - {tvK) g r ) (C8b) 
Ein(P n - , P n ) = r Ein 



1 



e{n)[^{{tvKf + \\Kf)g r -{txK)K 
+ L n (K -(tYK)g r )) . 



(C8c) 



APPENDIX D: Spherical symmetry 

We recall that the isometry group, Isom(A / l,g), of a 
spacetime (A4 , g) is the subgroup of the diffeomorphism 
group of A4, Diff(jM), which leaves the metric g invari- 
ant: Isom(7W,g) := {</> G Diff(M) | 4>*g = g}. 

Definition3 (Spherical symmetry). A four-dimensional 
Lorentzian manifold (Ai,g) is said to be spherically sym- 
metric if its isometry group, Isom(.M, </), contains a sub- 
group G with the following two properties: (i) G is iso- 
morphic to 50(3) and (ii) each orbit of G is spacelike 
and two-dimensional (up to some closed proper subset of 
fixed points). A tensor field T on a spherically symmet- 
ric spacetime is said to be spherically symmetric if it is 
invariant under G, hence if <f>*T = T for all <j> € G. 

From this definition it follows (excluding the case 
where the orbits of G are diffeomorphic to the 
two-dimensional real-projective space) that a four- 
dimensional spherically symmetric Lorentzian manifold 
(M,g) can, at least locally, be expressed as a warped 
product Ai = B Xr S 2 between a two-dimensional 
Lorentzian manifold (B, gg), called the 'base', and the 
stan dard unit two-sph ere (S 2 ,g.q^) , called the 'fiber' 



(see (|Straumannl . 12004 ) and (|Q'Neilll . [l983] )V This means 
that, at least locally, the manifold is a product 



M =B x S 2 
and the metric is given by 

9 = TT*(g B ) - (i?o7r) 2 cr*(gr S 2) 



(Dl) 



(D2) 



Here, tt and a are the projections of Bx S 2 onto B and S 2 , 
respectively, and 7r*,cr* their pull-backs. The warping 
function R is nothing but the areal radius, since, for a 
point p G B, the area of the fiber p x S 2 is just AirR(p) 2 . 

In this contest, a vector field X on M. at some point 
(p, q) £ B x S 2 has then a unique decomposition X = 



34 



tang X + taiig2 X in a component tangent to the 'leaves' 
B x q = a^ 1 (q) and a component tangent to the 'fibers' 
p x S 2 = 7r _1 (p). Arbitrary tensor fields on 2? and on S 2 
can be lifted to tensor fields on M. in the standard way. 
For covariant tensor fields (and hence, in particular, for 
functions) this is achieved via the pull-pack of the respec- 
tive projection: as an example, just look at (ID2j) . For 
contravariant tensor fields it suffices to consider the spe- 
cial case of vector fields. Let, for instance, X be a vector 
in the tangent space of B at p. Then the lift X at (p, q) 
of X is defined as the unique vector in the tangent space 
of M at (p,q) with tt^(X) — X and a*(X) — 0. Since 
this assignment is smooth, one gets the lifting of a vec- 
tor field via the pointwise lifting just described. In this 
work, we will mainly omit lifts and projections and not 
explicitly distinguish between original and lifted quanti- 
ties. For example, when referring to a vector 'tangent to 
B' we refer to a vector in the tangent space of B or to 
the lift thereof in the tangent space of M. . 

If X is spherically symmetric, then the component tan- 
gent to the fibers must vanish: tang2 X = 0. Similarly, a 
spherically symmetric one-form on M. must necessar- 
ily be tangent to B (i.e. normal to S 2 ) and thus it can 
be written as 8 = ir*(6g), where 6g is a one-form on B. 
Finally, a spherical symmetric function is simply the lift 
of a function on B. 



1. Connection and curvature decomposition 

In the following we discuss relations which express the 
curvature of M. in terms of the warping function R and 
the curvatures of the base B and the fiber S 2 . We start 
out from the relations between the Levi-Civita connec- 
tion V of (M,g) and the Levi-Civita connections of the 

B S 2 

base and the fiber, denoted by V and V, respec- 
tively. These relations can be derived, f or example, by 
means of the Koszul formula (see e.g. ( O'Neilll . Il983l ) 



Proposition 7.35). Let in the following X, Y ', Z be vec- 
tor fields tangent to B and U, V, W tangent to S 2 . Sup- 
pressing lifts and projections, we have: 

V X Y= B V X Y (D3a) 

Vx V = V V X = R- 1 X(R)V (D3b) 

tan S 2 V V W = 52 V v W (D3c) 

tang V v W = -g(V, W)R~ l VR. (D3d) 

Note that, for a function / on B, the lift of the gradient is 
equal to the gradient of the lifted function, that is (sup- 
pressing the lifts): grad/ = B grad/. For brevity, we 
write just grad / for it. Take care that for the Hessian 
and the Laplacian this is in general not true (see (|D13[) 
and (|D14|) V Therefore we write explicitly the super- 
scripts i B' in B Hess/ and B A / to denote the Hessian 
and Laplacian of / on B, respectively, or the lifts thereof. 

By means of (|D3[) one can now compute the expres- 
sions for the Riemann tensor. As the sectional curvature 



of S 2 is obviously constant and equal to one, the Rie- 
mann tensor, the Ricci tensor, and the Ricci scalar of S 2 
are simply given by s Riem = g s 2 g s 2 , s Ric = g s 2 , 
and s Seal = 2, respectively. Here we made again use 
of the Kulkarni-Nomizu product (|A8j) . Moreover, since 
the basis manifold B is two-dimensional, one can express 
its curvature tensors in terms of the scalar curvature. 
The expression for the Riemann tensor, Ricci tensor, and 
Ricci scalar of a spherically symmetric Lorentzian mani- 
fold (|D1ID2|) are, respectively, 

B Scal 

Riem = — - — g B Qg B 



1 + (dR, dR) ) R 2 g S 2 R 2 g S 2 (D4) 



I 

R 2 

2Rg s 2Q B HessR : 



B Scal 2 „ 

Ric = q B Hess R 

2 y R 



l + (dR,dR}+R e AR)g s 2, (D5) 



and 



Seal = B Scal - — f 1 + (dR, dR}) — — B A R . (D6) 
R d \ ) R 

Hence, for the Einstein tensor we have the expression 
Ein=ri(l+(dR,dR)) + ! B Ai?j g B - | B Hessi? 



3 Seal 1 



R 



A R R g s 



(D7) 



and for the Weyl tensor, using (|A9[) and (|A10aJ) . the sim- 
ple expression 

Weyl = w (g B g B + g B R 2 g s 2 + R 2 g s 2 R 2 g s 2 ) . 

(D8a) 



Here we put 



w-= \ ( B Scal--^(l + (di?,di?)) + -| B A R J (D8b) 



iscal+- B A R . 
6 R 



(D8c) 



In the derivation of (|D8|) we made use of the formula 
hsGgs = (l/2)(tr gB h B )g B 0g B , valid 22 for any symmet- 
ric bilinear form h B on B, in order to express the only 
term involving B Hessi? in terms of the Laplacian of R. 
Note that from (|D8a[) it is immediate that the Weyl ten- 
sor of a spherically symmetric spacetime has only one 



22 To prove this just note that the only independent component of 
this formula is the (eo, ei, eo, ei) one, where {e M } is an adapted 
orthonormal basis of (A4, g) such that eo, e.\ are tangent to B and 
e.2 , B3 are tangent to S 2 . Then, the equality follows immediately 
using the definition iASl of the Kulkarni-Nomizu product. 
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independent component, as it must be the case due to it 
being of Petrov-type D. 

Comparing expression jD4| with the definition of sec- 
tional curvature (|A6[) one can immediately read off that 
the sectional curvature K of the plane tangential to the 
(two-dimensional) S'0(3)-orbits at a given point is 



K 



1 

R 2 



(1 + (dR, dR)) 



(D9) 



and hence the MS energy, defined as (minus one-half) 
K times the third power of the areal radius, is given 
by ([58]) . Using the MS energy ([58]) . we can write the 
Einstein tensor (ID7D as: 



(D10) 



'Seal 



1 

R 



AR\R 2 g s 



We conclude giving the decomposition for the diver- 
gence of a spherically symmetric vector field (that is a 
vector field X tangent to B) and for the Hessian and 
Laplacian of a spherically symmetric function (that is a 
function /onK). First, from (|D3a[) and (|D3b[) . we obtain 
the following decomposition for the covariant derivative 
of X (expressed (0,2)-tensor): 



VI= B V X X{R)Rg S 2 . 



(Dll) 



Note that the mixed term (B-S 2 ) vanishes — as it should 
due to spherical symmetry. Taking the trace of (|D1 1[) 
one obtains the following expression for the divergence: 

div X = divg X + -^X(R) = div B (R 2 X) . (D12) 

The decompositions for the Hessian and Laplacian of a 
function / on B follow from inserting X = grad / in the 
above formulae. One gets: 

Hess/= B Hess / — g B (grad /, grad R)Rg S 2 (D13) 

and 

A/= B A/ + 2g B (grad/,gradi?)/i?, (D14) 
respectively. 

2. Einstein equation in case of spherical symmetry 

A general spherically symmetric matter energy- 
momentum tensor has the form 



T = T B +pR 2 g S 2 



(D15) 



where p is the spherical part of the pressure. Hence, using 
the decomposition (|D10[) of the Einstein's tensor found 
in Appendix [Dl the Einstein equation takes the form 

l{jL+ B AR\ g B -| B Hess R = 8nT B (D16a) 



B Scal 1 BA n a 

A R — 8irp . 

2 R F 



(D16b) 



Using the trace of the first equation, 

i(**+g)-*r«k. 



(D17) 



to eliminate e A R, one can write (|D16[) in the equivalent 
form 



1 f E 



R \R^ a ' 
'Seal 2E 
~Y~ + R* 



'Hess i? J = -47r*T e * 
= 4 7 r(trT jB + 2p) . 



(D18a) 
(D18b) 



Here, and in the following, * denotes the H odge-duality 
map for (B,g B ) (for the definition, see e.g. ( Straumannl . 



2004)). In the first equation we used the identity *t* = 
t — tr(x) g, which is valid for any bilinear form t on B, 
where the first (second) star acts on the first (second) 
slot of r. 

Finally, the integrability condition divT = for the 
energy-momentum tensor (|D15|) reads 



diw B (R 2 T B )+pd(R 2 ) = 0. 



(D19) 



3. Misner-Sharp energy 



Let us now turn to the MS energy and its properties. 
We first show that it is the charge of a conserv ed current . 
The treatment presented here follows mainly (jHavwardl 
Il996f ). In a spherically s ymmetric spacet ime one defines 
the Kodama vector field (|Kodam a. 1980) as the (unique 
up to a sign) spherically symmetric vector field orthogo- 
nal to, and of the same norm as, the gradient of R; hence 
we put 



k := -kdR . 



(D20) 



With this sign choice k is future-pointing if the gradient 
of R is spacelike. The orthogonality between k and the 
gradient of R is simply expressed by 



k(R) = 0, 



(D21) 



which clearly means that the integral curves of k stay at 
constant areal radius. 

An immediate but important property of the Kodama 
vector field is that it is conserved: 



div k = . 



(D22) 



Indeed, using (|D12|) and (|D21|) . one has: divfc = 
div/? k = dk = — *d-k -kdR = 0. 

Now, a key point for the study of spherically symmet- 
ric spacetimes is the following equation relating the MS 
energy with the matter's energy-momentum tensor: 



dE = 4nR 2 -kj 



(D23) 



where j is the so-called Kodama current (tangent to the 
base manifold B) defined by 



T(fc,-). 



(D24) 
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Equation (|D23[) follows from Einstein's equation; more 
precisely, it is equivalent to its fi-part (that is Eq. (|D18a|l ) 
fed with Vi?. To see this, just compute the differential 
of (5EJ as follows: dE = (E/R)dR + (R/2)d(dR, dR) = 
(E/R)dR + R B Hess R-dR = R((E/R 2 )g B + B Hess R) • 
dR, where the dot denotes here the contraction of the 
last slot of the tensor on the left with the first slot 
of the tensor on the right of the dot. In the second 
step we use that d((dR, dR))(X) = X(g(VR,VR)) = 
2g B (V x Vi?, Vi?) = 2 B Hess R{VR ,X) for any X tan- 
gent to B. Then, inserting (|D18a[) and using that * 
is skew-adjoint on one-forms, one gets dE = —AirR 2 * 
Tg * -dR = 4irR 2 * Tg • *dR and hence, using the defi- 
nitions (|D20j) and (|D24j) together with the symmetry of 
T, one arrives at (|D23|) . 
From (|D23[) it is clear that 



j(E) = 0, 



(D25) 



which means that the vector field j is tangent to the 
curves in B (hypersurfaces in M) of constant MS energy. 
Moreover, (|D23p implies that j is also conserved 



divj = 0, 



(D26) 



where the divergence is here taken on the spacetime 
(M,g). To see this, just compute the divergence of j 
with (|D12p and using the Hodge-dual version of (|E)23|) : 
divj = R- 2 div B (R 2 j) = R- 2 S(R 2 D = (inR 2 )-^ * 
dE = 0. 



Following ( Havwarcj 1 1 996} ) , we can now show that the 
charges corresponding to the conserved currents j and 
fc are, respectively, the MS energy and the areal vol- 
ume. Let S be some spatial three-dimensional hypersur- 
face which, because of spherical symmetry, decomposes 
as £ = a x S* 2 , where a is some spatial curve in B. Recall 
that the charge related to a conserved current X is given 
by Qx(S) ■— f^i-xl^, where fi is the volume form on 
M, and, because of spherical symmetry, the latter de- 
composes as /it = hb /\R 2 (is 2 , where /xg and n$ 2 are the 
volume forms on B and on the unit two-sphere, respec- 
tively. After integration of the spherical part and since 
ijfJ-B = *j, using (|D23j) one gets 



<&(£)= / dE, 



(D27) 



which means that the charge of j is the MS energy This 
justifies the interpretation of the MS energy as a quantity 
associated to the 'interior' of the considered sphere of 
symmetry. In fact, due to (|E)26|1 . the charge does not 
depend how one choose the spatial slice to define the 
interior. Similarly, since ikfJ-s = *fc = dR, the charge to 
fc is simply J AirR 2 dR and hence 



Qik(S) 



(D28) 



which says that the charge of fc is the flat-space volume 
computed with the areal radius. 



Incidentally, the Kodama vector can be used to give 
an elegant proof of Birkhoff 's theorem, which states that 
spherically symmetric solutions of Einstein's equations 
are, in fact, static. Indeed, by direct computation one 
shows that in vacuum fc is Killing and, because of spher- 
ical symmetry, it is clearly also hypersurface orthogonal. 

Next we turn to the relation be tween the MS e nergy 
and the Hawking quasi-local mass (|Hawkinel . ll968l ). The 
latter is a quantity associated to a spatial two-sphere, S, 
in an arbitrary spacetime. It is defined by 



M H (S) 



Arez(S) 
16tt 



f 



f 

2^ 



Ms 



(D29) 



Here, 9 ± := tr S 2(Vi ± )/2 are, respectively, the expan- 
sions of the outgoing and ingoing future-pointing null 
vector fields I normal to S, the latter being partially 
normalized such that g(l + ,l~) — 1 (there remains the 
freedom to rescale I — > a^ 1 ^, where a is a positive 
real- valued function). In the special case of spherical 
symmetry we take S to be an orbit of the rotation group. 
Then we clearly have Area(S') = AttR 2 . It is also obvious 
that the metric of the base B, evaluated on S, can simply 
be written in the form 



z + <g)Z +1 <s>r 



(D30) 



Now, for V tangent to S, (|D3b)) gives Vyi 1 = 
R- 1 l ± (R)V so that 6 ± = i?- 1 Z ± (i?). Hence we have: 

2 6+0- = 2R~ 2 dR(l + )dR(r) = (dR, dR) /R 2 , (D3I) 



where we used (|D30I) , or rather its contravariant version, 
in the last step. Equation (|55|) now establishes the equal- 
ity between the MS energy at p and the Hawking quasi- 
local mass of S, where p is any point on S: 



E(p) = M U {S) 



(D32) 



As is the case for the Hawking quasi-local mass, the 
MS energy can be naturally decomposed into a Ricci and 
a Weyl part: 



where 



E — E R + E v 



E R :— —^R 3 K R , 



Ew 



(D33) 

(D34a) 
(D34b) 



Here K R and K w denote, respectively, the Ricci and the 
Weyl parts of the sectional curvature of the plane tan- 
gential to the S l O(3)-orbits. These are obtained inserting 
the decomposition of the Riemann tensor (|A9|) in the def- 
inition of the sectional curvature (|A6|) . The Ricci part of 
the MS energy is determined by the local matter distri- 
bution via Einstein's equation: Using expressions (|D15|) 
and (|D2p for an arbitrary spherically symmetric energy- 
momentum tensor and, respectively, metric in (|A12|) one 
gets 



E„ 



-in 



R 3 (trT B +p). 



(D35) 
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For the Weyl part of the MS energy we have, in view 
of (jD8a|) . that 

E m = -±R 3 w, (D36) 

where w is given by (|D8b[) or (|D8c[) . Hence, in particular, 
the Weyl tensor vanishes iff E w does. Since the square 
of the Weyl tens or is (Weyl, Weyl) = W aM& W a ^ & = 
12w 2 , with (|D36[) we obtain the nice expression 

(Weyl, Weyl) =48^. (D37) 

From this one sees that, in a spherically-symmetric space- 
time, the non-vanishing of E m (that is the non-vanishing 
of the Weyl tensor) for R — ► implies a curvature singu- 
larity at R = 0. 

To gain a better physical understanding of the Weyl 
part of the MS energy we take a look at the equation 
of geodesic deviation (sometime called Jacobi equation). 
Let u be a geodesic observer field and s the spatial 
{g{u, s) = 0) separation vector between two nearby inte- 
gral curv es of u. Then the e quation of geodesic deviation 
(see e.g. (|Straumannl . 120041 )) is V„ V u s = R(u, s)u = 
B u (s) + C u (s). In the last step we decomposed the en- 
domorphism on the r.h.s. of the geodesic deviation equa- 
tion in its Ricci- and Weyl-part, denoted here by B u and 
C u , respectively. For an arbitrary spherically symmetric 
spacetime the latter is given by (see (|D8|) and (JD36J)) 

C U = —j^fPu - ~]^ P U > (D38) 

where P° and P£ 2 are, respectively, the B and S 2 
parts of the projector P u . Recall that P u projects onto 
the subspace of the tangential space orthogonal to u 
(see (|125bjl ). Equation (|D38|) is exactly the same expres- 
sion one gets in Newtonian gravity — provided one identi- 
fies E w with the mass of the central object. The spatial 
endomorphism C u just describes the familiar volume- 
preserving tidal deformation which produces an expan- 
sion in radial direction and a contraction in the orthogo- 
nal directions tangential to the S'0(3)-orbits. 

Concerning the Ricci part B u , in the case where u 
is the velocity field of dust 23 , making use of Einstein's 
equation (see (|A12jl ) we have: 

47T 

B u = - — qP u . (D39) 

This just says that the local effect of matter (here given 
by dust) is an isotropic contraction. 

4. Spherically symmetric perfect fluids 

We specialize now to a perfect fluid, which is described 
by a four-velocity vector field u, density g, and pressure 



Recall that dust particles moves along geodesies by the Euler 
equation. 



p. In case of spherical symmetry u is tangent to the basis 
manifold and the matter energy-momentum tensor (|70p 
decomposes as 

T = gu®u+p(u®u- g B ) +pR 2 g s 2 , (D40) 

from which one can read off the part tangent to B: 

T B = gu(g>u+p(u(g>u~ g B ) . (D41) 

Usually, the description is to be completed with the spec- 
ification of an equation of state. We will not assume any 
equation of state yet, since in some cases (e.g. McVittie 
spacetime) this happens to be determined by Einstein's 
equation. 

Inserting (|D41[) in (|D35[) we get for the Ricci part of 
the MS energy the simple expression 

47T 

E r =—R 3 q. (D42) 

Also the expression (|D23[) for the differential of the MS 
energy simplifies in case of a perfect fluid. Using (|D41|) 
the Kodama current (as one-form) becomes 

1= {g + p)g{k,u)u- pk. (D43) 

It is useful to introduce an adapted orthonormal basis 
{it, e} tangent to the basis manifold, where u is the ve- 
locity vector field of the fluid and e is chosen to point 
in direction of increasing areal radius. Because of our 
choice of orientation we have e = -ku (the volume form 
on B is simply fig = uAe). Using this expression for 
e and the definition of the Kodama vector field (|D20j) 
we have g(k,u) = (fe,u) = {*dR,u) = — {dR, *u) = 
— (dR,e) = —dR(e) and hence, the Hodge star of the 
Kodama current becomes 

★ j = -p dR(u)u - g dR(e)e , (D44) 

which, inserted in (|D23|) . gives the following expression 
for the differential of the MS energy for a perfect fluid: 

dE = -4vri? 2 (p dR(u)u + g dR{e)ej . (D45) 

Hence, the variation of the MS energy along it and e is, 
respectively: 

dE(u) = -4irR 2 pdR(u) , (D46a) 
dE(e) = +4ttR 2 g dR{e) . (D46b) 

These expressions have a good physical interpretation: 
Since the matter moves along u, (|D46a[) expresses the 
fact that the energy can only increase (decrease) if the 
motion along u does (releases) work against (with) the 
action of the pressure. Equation (|D46bj) expresses the 
almost obvious increase (decrease) of gravitational mass 
with increase (decrease) of volume in the rest system of 
the matter. We said 'almost' because A-KR 2 dR(e) is not 
quite the increment of proper volume. The difference ac- 
counts for the fact that kinetic and gravitational binding 
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energy are themselves gravitationally active. To see that 
this is indeed what (|D46b[) implies, let p be some point in 
spacetime and S p the two-sphere of spherical symmetry 
through p. Assume S p to have a regular interior, that 
is, that S p bounds a 3-ball B p in the hypersurface E or- 
thogonal to u. Except for the origin of B p , we can write 
B p — a x S 2 , where a is a spacelike curve in B orthogonal 
to u, going from the center of symmetry to w(p). Using 
the expression E = (R/2)(l + (dR(u )) 2 - ( dR(e)) 2 ) for 
the MS energy to eliminate dR(e) in (|D46b[) . integrating 
the latter over a, and re-expressing the result as a volume 
integral, one gets: 



and the shear scalar 



m = I 

Jb, 



LI {l + (dR(u)) 2 -^ 



1/2 



MS 



(D47) 



One sees that the MS energy contains the contribution 
from the proper mass contained in the ball B pi 



M{p) = I g fXY, i 



(D48) 



as well as c ontributions from the 'kinetic' and 'poten- 
tial' energy (jHavwardL Il996t iMisner and Sharp! . 1 1964 ). 
In a Newtonian approximation, that is for small 'veloc- 
ity' dR(u) and weak field (small E/R) one can expand 
the square root in (|D47j) and gets, in leading order: 



E(p) 



h(dR(u)) 2 



R J 



Ms 



(D49) 



In this approximation the MS energy is therefore just 
the sum of the proper mass and the Newtonian kinetic 
and potential energies contained in the ball B p . This 
provides a sound justification for the interpretation of 
the MS energy as the active gravitational energy. 

At this point we can compute also the differentials of 
the two parts (|D34[) of the MS energy separately. The 
differential of the Ricci part follows directly from (|D42[) : 



dE R = AnR 2 (g dR+\R dg) 



(D50) 



and the differential of the Weyl part is just the difference 
of this with (|D45]) : 



AnR 2 1(q + p)dR(u) u+\Rdg) . (D51) 



dE v 



Its components in the directions u and e are then 

dE w (u) = -AnR 2 (g + p)dR(u) - ^-R 3 dg(u) , (D52a) 
<2£ w (e) = -^R 3 dg(e). (D52b) 

It is now instructive to express the variation along u 
of the Ricci and Weyl parts of the MS energy in terms of 
the kinematical properties of the fluid velocity u. Recall 
that, because of spherical symmetry, the rotation tensor 
vanishes identically and the shear tensor has only one 
independent component. The kinematical quantities re- 
duces thus to two scalars: the expansion 



dR{u) 
R 



(D54) 



The shear tensor is then given by the trace- free endomor- 
phism <t = a(Q S 2 — 2 Q e ), where Q S 2 and Q e denote, re- 
spectively, the projections onto the two-dimensional sub- 
space of T(M.) tangential to the two-sphere and onto 
the one-dimensional space parallel to e (for the latter 
see (jClaP ). We recall that the divergence-freeness of the 
energy-momentum tensor (|D40[) is equivalent to 



(g + p)6 = -dg{u) 
(g + p)b = -dp(e) 



(D55a) 
(D55b) 



where b := — <jf(V„it, e) is the acceleration (scalar) of u 
in positive radial direction (the minus sign in the lat- 
ter formula is because the metric is negative definite in 
spatial directions). 

Now, using (|D55ap and (|D54[) we get: 



47T 

dE R (u) = —R 3 (3ga-p9) 

47T 

dE w (u) = -—R 3 (g + p )3a 



(D56a) 
(D56b) 



9 := div u 



(D53) 



With the equations just derived we can now say when 
the MS energy, and its Ricci and Weyl parts, are tem- 
porally or spatially constant. Here, by temporally (spa- 
tially) constant we mean that the variation in direction 
of u (e) vanishes. We collect the results in the following 

Theorem 4. Consider a spherically symmetric fluid with 
g + p and restrict to the region where dR is spacelike. 
Then for the MS energy E and its Ricci and Weyl parts 
E R and E w the following statements hold true: 

(i) E is temporally constant iff p = or dR{u) = 0; 
(ii) E is spatially constant iff g = 0; 
(Hi) E w is temporally constant iff a = 0; 
(iv) E w is spatially constant iff g is spatially constant; 
(v) E R is temporally (spatially) constant iff R? g is tem- 
porally (spatially) constant. 

The proof is a straightforward application of the formulae 
just derived above. Note that the assumption g + p =/= 
is needed only for (Hi). The assumption that dR is 
spacelike is needed only for (ii): If dR is spacelike, then 
for any spacelike spherically symmetric vector e (hence 
tangent to the basis manifold B) it holds dR(e) ^ 0, 
since in a two-dimensional Lorentzian manifold any two 
spacelike vectors are linearly dependent. 
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